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OPTIMAL LIMIT STATE DESIGN OF 
REINFORCED CONCRETE 
SHAFT SUPPORTED CONICAL WATER TANKS 

(A thesis submitted in partial fulfilment of the requirements 
for the degree of Doctor of Philosophy by M.L. Kalwar to the 
Department of Civil Engineering, Indian Institute of Technology 
Kanpur,, India) 


Water towers are the land marks of any modem day 
habitat. Unlike in the past, the present day society expects, 
rather demands the modem structures to be shaped in such a way 
that whole environment is enhanced by their very presence. In 
a developing country like India wherein the resources are 
limited, a structural designer has to itake a compromise between 
the aesthetic and economic aspects. A conscientious designer 
has to proportion the structure in such a fashion that struc- 
tural efficiency and economy are acdiieved simultaneously. 

Hence the role of a present day structural designer is all the 
more exacting as he has to design a structure which satisfies 
the requirements of safety, functional aspects, economy and 
aesthetics. The choice of a structural material depends upon 
its initial and the maintenance costs. Hence reinforced 
concrete water towers are preferred over others. A water tower 
comprises of a tank or a reservoir, a supporting structure and 
tlte foundation. The present study addresses itself to the 
superstructure only viz, the tank and its supporting structure 
The water towers considered in the present study are 
multishell structures consisting of a conical tank supported 



on a cylindrical ^aft with a spherical top. These towers are 
classified as Type I, II and III depending upon the capacity 
and the internal geometric arrangement of the elements in the 
reservoir portion, A multi^ell structxire may be composed of 
simple shell and circular ring elements. The elements consid- 
ered herein are axisymmetric in nature. 

Traditionally the reinforced concrete structures are 
designed either by working stress method based an elastic theory 
or by \iltimate load theory. Both of these methods are found to 
be inadequate to provide a rational design. It is also realised 
that the loads and the material properties are after all not 
deterministic. The foregoing considerations set the stage for 
the development of a much needed, more rational and sovind 
design philosophy. The limit state design concept is a conse- 
quence of this. With the evolution of this new concept of 
design both the elastic and plastic methods of analysis together 
are brought into limelight. The codes of practice around the 
Globe are being revised to incorporate the limit state theory, 
and Indian Standards Institution is no exception. 

In the present study, the optimal limit state design 
of reinforced concrete, shaft supported, conical water tanks is 
investigated using the limit state method for Indian environ*- 
mental conditions. In India many water supply schemes are 
coming up and obviously a large number of water towers will 
have to be constructed to meet the requirement. Hence optimal 
design of such structures will be very meaningful. 

The response of a structure to service loads is 
predicted by elastic analysis and the resistance to ultimate 



loads obtained through limit or plastic analysis. The design 
forces in the reservoir portion are found to be governed by 
a xi symmetric loading combination of dead load and water loading 
xinder full tank condition. For this loading combination elastic 
analysis of the tower is carried out by tasing force method of 
analysis, which is specially suited to a xi symmetrically loaded 
^ells of revolution. However lateral forces due to wind and 
seismic effect under full or empty tank condition are found to 
be critical in case of supporting shaft. Hence for simplifi- 
cation the shaft is idealised as a cantilever for purposes of 
analysis. For computing the seismic forces, response spectrum 
method has been used. 

For the proper application of the limit state theory, 
the resistance of the structure to ultimate loads is to be 
ascertained. Therefore, various collapse mechanisms of rein- 
forced concrete conical shell subjected to a hydrostatic 
pressure, adopting the upper bouxxi approach, are investigated. 
All the possible mechanisms are considered to arrive at an 
exact value of collapse pressure. In Type III water towers, 
in addition to conical reservoir, the collapse mechanisms of 
cylindrical partition are also investigated to arrive at the 
collapse vgater pressure. The resistance of the supporting 
shaft to ultimate loads is determined by carrying out its 
limit analysis. 

All the strength and serviciability constraints of 
the relevant limit states are incorporated in the nonlinear 
programming model of the formulated optimum design problem. 
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•Thie shaft diameter# ratio of top diameter to bottom diameter 
•of the conical reservoir# thicknesses of various elements and 
■HJi-.e percentages of reinforcement in these elements# with 
Kiuacticable logical distribution of reinforcement, are taken 
SIS d-esign variables. The overall cost of the superstinacture 
cof tie tower# which inclxxies the costs of material^ formwork 
sncd labour# is c±iosen as the objective function for the 
copftimal design. The optimization study of these towers yields 
nott only the optimal value of the design variables but also 
tthse optimal configuration of the tower. Such an investigation 
Inaaslieen carried out for the first time. Furthermore, the 
seras-ltivity of the optimal design with respect to tiie escala- 
ticaa of tinit costs of concrete, steel or formwork has been 
stuudded in the present work. 

Parametric studies have been conducted for the most 
csotEWtanly used reservoir capacities and staging heights. In 
coicBe*: to obtain designs for various wind zones of India, the 
boasic wind pressure is also treated as a parameter in the study 
oof optimal design of these water towers. The suitability of 
tlliesse designs, obtained for various wind zcxies, are establ- 
l-shceSl for the different seismic zones. 

In order to save the ccanputaticmal time the optimi- 
asatcioii is carried out using an indirect approach to arrive at 
L-hes optimal configuration. This approach has not only resul- 
Ue* in ccnsiderable saving of ccanputational time but also 
cpdws designs with optimal design variables and the optimal 
r.'attio of top diameter to the bottom diameter of the conical 
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reservoir for different ^aft diameters, in addition to the 
one which corresponds to the optimal configuration. 

For the benefit and easy reference of the field 
engineers, designs are presented in tabxiLar form from which 
the optimal values of design variables may be picked up for a 
particular tank capacity, staging hei^t, wind and seismic 
zones. Thus the designer will be able to arrive at an optimal 
design which conforms to the relevant codes of practice and 
without recourse to a powerful computing device. The compu- 
tations have been carried out with the standardisation of 
form in view. Thxxs, the optimal solutions presented are the 
realistic solutions, embodying the field considerations, 

Althoui^ the present study is confined to Indian 
environment, the results can be judiciously utilized for other 
environments as well. In essence the thesis contains optimal 
limit state designs of reinforced concrete shaft supported 
conical v«iter tanks, which are useful to a professional 
engineer. 



CHAPTER 1 


INTRODUCTION 


1.1 General 

The basic aim of a stnjctural engineer is to design 
such structviral systems which satisfy safety, fxanctional # 
aesthetic, and economic requirements in the most efficient 
manner r in order to achieve the best results, one has to adopt 
such a methodology, for analysis as well as design, which is 
capable of predicting the realistic behaviour of a structure 
when subjected to variovis design loads and gives satisfactory 
design which remains safe and serviceable throughout its 
design life. However, if a structural theory is req\iired to 
take into account every possible variable involved, it will 
become too complicated for practical use. So the general 
practice is to make simplifying assumptions that yield consis- 
tent and sufficiently accurate results. Experience, experi- 
ments and baaig understanding often are required to determine 
whether a giv^ theory or method is applicable to a particular 
situation. Neither 'working stress design method' nor ' ulti- 
ma te load design theory' has been fovind to be adequate for 
design of reinforced concrete structures. Moreover, to take 
into account the variaticais that are likely to occur in the 
loads to which the structures are subiected and the varia- 
tions in the strength of material of which they are comprised, 
a probabilistic approach to safety has to be adopted. 



No doubt/ ‘limit state design philosophy' eliminates 
the limitations of working stress and ultimate load design 
methods and is capable of giving satisfactory designs. But, 
in order to ensure economy and efficiency in such designs, 
optimization techniques have to be incorporated in the design 
process. 

Recent advances in computer-aided analysis and 
design of concrete structures have made it possible to analyse 
and design complex concrete structures xander the ejtpected 
static and dynamic loadings. The computer is an important 
aigineering tool allowing considerable automation in the 
standardised design of conventional concrete structures. 
Optimization becomes more meaningful when used for such 
standardised designs, 

Ihe availability of a very good computer facility at 
the Indian Institute of Technology , Kanpur and the foregoing 
considerations have motivated the present study dealing with 
optimal limit state design of reinforced concrete water towers. 

1,2 The Present Study 

In a developing cotintry like India one of the most 
important basic needs is to provide adequate protected water 
supply systems to industrial and residential complexes in 
cities, towns and villages. Water towers are among the major 
constructicais of most water supply systems. Many water supply 
schemes are coming up throughout the country and obviously a 
large nuinber of water towers will have to be constructed to 
fulfil the need. In the recent past many reinforced concrete 





water tcwers/ particularly with cylirvirical shell type support, 
have collapsed in India. This has resulted in some of the 
organisations losing faith in the construction of shaft over- 
head tanks. On investigation of some of these failures, 
it was observed that a number of reasons are responsible for 
this sorry state of affairs. Ihe principal reasons identified 
are , 

(i.) lack of uunder standing of the behaviour of the stmc- 
tural system; 

(ii) the ignorance of the tolerance limits on the geometric 
parameters; and 

(iii) the use of substandard formwork in the construction 
stage. 

In view of these facts a deeper insight into the structural 
behavioxar and design of such structures is sought for. 

Moreover, not much of work has been done on optimization of 
water towers using either elastic or plastic models. 

In the light of the foregoing discussions, the obvious 
reasons for undertaking the present study may be summarized 
as: 

(i) To analyse and design reinforced concrete water 

towers using modem numerical tools which allow a 
useful insi<^t into the structural behaviour and 
design aspects of such structures under the expected 
static and dynamic loadings. As the number of such 
structures to be constjructed is very large, any 
small saving that can be made in the design will get 


added up. Hence optimization technique is incorpo- 
rated in the design process. 

(ii) To provide, for the benefit of the field engineers, 
optimal standardised designs of reinforced concrete 
water towers with various capacities and staging 
heights for different zones of India based on wind 
and seismic forces. This will make the designers 
free from substantial tedious work. Besides, stan- 
dard forms may be designed for repeated use. 

(iii) To instil interest among the academic community that 
such optimal studies usually provide, specially in 
areas where not much work had been done. 

1,2.1 Choice of the shape for the water towers 

Water towers are among those structures which come 
under the category of prominent features of a particular 
locality. Besides economy the aesthetic value of such struc- 
tures is also of considerable importance and hence a proper 
consideration should be given to this aspect while choosing 
a particular shape for such structures. Keeping this fact in 
view, a conical reservoir with a cylindrical support has been 
chosen, in the present study, as the design configuration for 
the water towers. Such a configuration, undoubtedly, is 
aesthetically pleasing and has better architectural expre- 
ssiveness than Intze tank or cylindrical tank supported 
either can columns or :^:iaft. Moreover, this configuration 
does not pose any extra difficulty in construction. 



1,2.2 Capacities, staging heights and the types of tanks 
considered 

The scope of the present study has been kept limited 
to the designs of reservoirs of capacities ranging from 100 
kl to 1000 kl on stagings of 15 m, 20 m and 25 m heights as 
these will be among the most commcaily used capacities and 
heights. Reservoir capacities considered in the present stxidy 
are 100,;-200, 300 , 500 , 750 and 1000 kl . Depending upon the 
capacity and other requirements the internal geometric 
details of the reservoir are fixed without any change in tine 
overall configuration. In the present study, based on the 
internal details of reservoirs, these water towers have been 
classified into three categories as Type I, Type II and 
Type III, The Figure 1,1 shows the details of these towers. 

Type I — 

This Category ^own in Figure 1.1(a) is intended for 
the reservoirs of capacities 100 to 300 kl . Reservoirs upto 
300 kl capacity require supporting shaft of small diameter, 
Hence^ in this case, diameter of internal i^aft is kept same. 
Ihls will result in easy constructional details without losing 
much of space of the reservoir. 

Type II — 

For reservoirs of capacities 5(X) to 1000 kl, towers 
of the type shown in Figure l,l(b) are usually preferred. 

Xs the diameter of supporting shaft required for these capa- 
cities will be more, reservoirs are provided with spherical 
bottc»T» and internal ^laft having smaller diameter as compared 
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to that of the supporting shaft. This geometric arrangoiient 
will provide more space for storage of water therdDy leading 
to overall eccsnomy. 

Type III — 

Sometimes a reservoir is required to have a partition 
for various reasons. The Figure 1.1(c) shows a reservoir of 
this type in which a cylindrical partition divides the reser- 
voir in two parts. Water towers of this type are also consi- 
dered for 500 to 1000 kl capacities. 

ThuSj. in brief/ the present study deals with the 
optimal limit state design of reinforced concrete water towers 
consisting of conical reservoir supported on cylindrical shaft 
for Indian aivironmental conditions, 

1.3 Limit State Design 

For the design of reinforced concrete structures# in 
general/ working stress method ax>d ultimate load theory have 
been extensively used. Both of these methods are determin- 
istic in nature. In the elastic or working stress apEaroach# 
the structure is assumed to act as a linear elastic body under 
the action of specified loads. The members of the structures 
are so proportioned that the maximum stresses in concrete and 
steel/ produced in the most critical section of the structure# 
for the worst corabinatico of loads# do not exceed the permi- 
ssible xrailues. These permissible values are obtained as a 
fraction of the specified strengths of ccncrete and steel. 

Many structures have bean successfully designed by this method 
The objections to this method are not that it is unsafe# but 



that it is idiysicaliy unreal and leads unnecessarily to 
conservative designs with varying degrees of safety against 
collapse. The maximum stress may occur at only one point in 
the whole str\x:ture which means that safety is being assessed 
by looking at a very local effect. It is also recognised 
that reinforced concrete deviates from elastic b^aviour 
appreciably ev«i at lower loads and radically at higher loads. 
Hence elastic analysis will be much less valid in case of 
reinforced concrete structures. 

Foregoing limitations of the elastic analysis led, 
in the post-war years, to the development of the earlier 
method of ultimate load theoiry which had been abandoned with 
the advent of the theory of elasticity. According to the 
xiltimate load theory, also termed as plastic or limit anal- 
ysis, failure does not occxir because of the stress reaching 
the limiting valxie at one point. On the other hand, it takes 
into accovmt the inelastic behaviour of the material and 
redistribution of stresses. Experimental results have con- 
firmed that the mode of collapse arxi load producing it are 
remarkably well predicted by this method. The margin of 
safety against collapse, called the load factor, can be 
Chosen as per the requirement. Another contribution of 
this method is that it has diverted the attention frcrni the 
microscopic stress concept to the macroscopic overall beha- 
viour of the structure. Althoii^ this is an excellent model 
for finding the straigth of a structure, it has a serious 
limitationi that it cannot predict the response of the system 
to service loads. It is observed that some of the most 



cotTinon causes of damage in reinforced concrete structures are 
excessive deflection, cracking etc, which make the structure 
mserviceable . Hence safety against these damages, deP!\and 
due considerations in the design. 

Thus, it is recognised that elastic theory is better 
suited to predict the response of a structure under service 
loads only arxi the ultimate strength can be determined only 
through limit analysis, \s pointed out earlier, both these 
methods are deterministic in concept. The randomness involved 
in the basic design parameters has rendered the deterministic 
approach obsolete. The foregoing considerations set the 
stage for the development of a much needed, more rational and 
sotand design philosojdmy. The limit state theory is a conse- 
quence of this and has been adopted for the design of struc- 
tures considered in the present study. 

1.3,1 Definition of limit states 

A structure, or a part of a structure, is rendered 
unfit for use when it reaches a particular state, called a 
‘limit state' in which it ceases to fulfil the function, or 
to Satisfy the condition, for which it was designed (TC98, 
1973), Limit states can also be defined as; A condition where 
a structure or component ceases to fulfil its intended func- 
tion in scxnre way, or becomes unfit for its intended purpose. 

Depending upon the safety, serviceability and other 
requiren^nts , limit states may be grouped into three major 
categories. 



(i) Limit states of strength: 

Those corresponding to maximum load carrying capacity 
and may be governed by the rupture of critical sections, or 
excessive plastic deformation leading to collapse, or instab- 
ility etc. 

(ii) Limit states of serviceability: 

Those corresponding to excessive deformations, wide 
cracking, excessive vibration etc. making the structure unser- 
viceable. 

(iii) Other limit states: 

Those corresponding to fatigue, lightning, fire, 
explosion etc. 

1.3.2 Safety formats for limit state design 

Various national and international cotrenittees, working 
with the aim to provide a rational and sound design philosophy, 
have recommended a semi- probabilistic limit state method. 

Thus, probabilistic concepts, for the first time, are expli- 
citly embedded in the design. It is assumed that the statis- 
tical distribution of loads and material strengths are avail- 
able, The ’characteristic values' for both loads and material 
strengths are then selected with acceptable probabilities 
against higher load or lower strength respectively. The 
practice of using a single global factor of safety has been 
replaced by probability-based partial safety factors, one for 
loads and the other for material strengths. The 'design 
values' are derived from the characteristic •'^lues through 



the use of partial safety factors. The partial safety factors 
for loads vary with the degree of seriousness of the parti- 
cular limit state being reached^ probability of two or more 
loads occurring together^ and the reliability of the struc- 
tural theory being used. Whereas, the partial safety factors 
for material strengths take into accotint the difference between 
the strengths of controlled specimens and the material in the 
prototype structure, with these modified loads and material 
strengths, safety of the structure is investigated for the 
relevant limit states. 

1.4 Optimi za t ion 

1.4.1 Structural optimization 

Any design philosophy demands that the components of 
the structure be so designed that performance requirements 
are satisfied. Such requirements are termed as behaviour 
constraints in optimization literature. Construction or 
manufacturing processes may also impose certain restrictions 
on the design which are termed as side ccnstraints. The 
usual method of indirect design, which does not make use of 
any optimization technique, is a trial and error process. 

An initial design is intuitively selected and the analysis is 
carried out to check its safety. If any constraint is vio- 
lated, or margin of safety is too high the design variables 
are suitably altered and the analysis is carried out again. 

This procedure is reputed till a satisfactory feasible 
design is obtained, kiwang the many possible feasible designs 


of a structxire, scxne designs are better than others and it is 
very difficult to choose the best one. On the other hand, 
the optimum seeking methods, also known as optimization tech- 
niques, start with an initial design and proceed systemati- 
cally to find the best design based on certain criterion, 
such as minimum weight/cost or maximum stiffness etc. More- 
over, the optimum design approach provides a rational basis 
for autcmatic design methods using a digital computer. 

1.4.2 Optimization techniques 

Basically there are two approaches to optimization, 
analytical and numerical. Usually analytical methods like 
variational calculus, Langrange multipliers etc. are employed 
to solve relatively simple optimization problems, with the 
advent of digital computers, nxamerical approach of optimiza- 
tion has become popular. These nxamerical approaches are 
available in the literature of mathematical programning and 
are capable of dealing with a variety of problems in the field 
of optimization. Linear, nonlinear, geometric, dynamic, 
quadratic and integer programming techniques fall under the 
category of mathematical programming. 

Linear programming deals with a particular class of 
optimization problems in which objective fxincticxi as well as 
constraints are linear functions of design variables. Most 
of the structxJTal engineering design problems do not come in 
this category. Both objective function as well as constraints 
are nonlinear in nature in most cases of structural design 
problems. Hence the most suited and widely used technique. 



for such problams^ is the nonlinear progran\ming one, which 
has been chosen for the current study as well. 

The nonlinear programming probl«ns can broadly be 
classified as 
(i) Unconstrained, or 
( ii ) con s tra ined , 

1.4.3 Unconstrained nonlinear programming 

the name implies, unconstrained nonlinear progra- 
mming deals with the optimization of an objective function 
without any constraint. This approach can be applied to 
constrained optimization problems also, after transforming 
the original problem to an equivalent unconstrained problem. 
Thus, mconstrained nonlinear programming has become one of 
the most important mathematical tools for solving both tancons 
trained and constrained optimization problems. 

Optimization is concerned with minimization or maxi- 
mization of an objective function. Since the maximvan of a 
function corresponds to the minimum of the negative of the 
same function, optimization can be stated as minimization of 
a function without any loss of generality . Thus , the uncons- 
trained optimization problem can be stated as finding the 
vector X of the design variables which minimizes the 
oblective function P(x) vhere there are no restrictions an. 
the choice of X. Xll the iinconstrained minimization methods 
are iterative in nature and hence they start frbm an initial 
trial solution ana proceed toxards the mlnlrtium point in a 
sequential manner. Ihus, if x^^ is the approxiAsltion to the 
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minimum in the (i-l)th iteration, the improved approximation 
in the ith iteration is found from the relation 



-f a* 


( 1 . 1 ) 


where 


= vector of design variables at the new point, 
= vector of design variables at the previous 
point , 

= a suitable direction of move along which the 
value of objective function decreases, and 
®i ^ step length which minimizes F(X) in the 

directicMi 


All the constrained minimization methods differ from 
one another c«ily in their approach to find and a*. These 
methods can be classified into two broad categories as direct 
search or descent methods. 

The direct search methods require only objective 
function evaluaticxis and do not use the partial derivatives 
of the functicffi in finding the minimum and hence are often 
called the nongradient methods. Random search method, 
univariate method, pattern search methods, simplex method 
are among the nongradient methods. These methods are, in 
general, less efficient than the descent methods. The 
descent methods are also Imown as gradient methods because 
they make use of derivatives of the objective function in 
addition to objjective function evaluation in finding the 
minimum, steepest descent method, FI etcher -Reeves method. 
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Nevsrton*s method^ Davidon-Pl«tcher--Powell method are among the 
gradient methods, 

Ihe various methods used to find the minimizing step 
length can be classified as elimination methods and inter- 
polation methods. Fibonacci method and golden section method 
among the elimination methods, and quadratic interpolation and 
cubic interpolation method among the interpolation methods are 
more contnonly used depending upon the mathematical nature of 
the optimization problem at hand. 

1.4,4 Constrained nonlinear prograitming 

Most of the engineering designs come under the cate- 
gory of constrained optimization problems. A constrained 
optimization problen can be stated in the standard form as: 

Find vector X of the design variables which minimizes 
the objective function F(X) subject to m constraints 

gj(X) < 0 f j “ 1# 2, ... m 

There are many techniques available for the solution 
of a constrained nonlinear programming problen. These methods 
can be classified into two broad categories, namely, the 
direct or indirect method. The methods of feasible directions 
and constraint approximation methods are among the direct 
methods. Indirect methods involve the transformation of 
ccaistrained problem into an equivalent unconstrained minimi- 
zation problem. The transformed problem is then handled by 
any of the methods of unconstrained minimization as a 

d ■ 
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sequence of mconstrained minimization problems. This is 
achieved either by transformation of variables or by using 
penalty function methods. 

1.4,5 Choice of optimization technique 

The choice for a particular method of optimization , 
in nonlinear progratiroing, depends on the size and the nature 
of the problem. The efficiency and the stability of these 
methods are also very much dependent on the problem. 

For the medium-size problem (10 to 50 variables) 
Davidon-Fletcher-Powell method, also known as variable metric 
method, has been found to be stable and efficient. Xs such, 
this method has beoa chosen for the current study to find out 
the search direction S. Since this method makes use of 
gradient of the objective function, it has been coupled with 
the cubic interpolation method, vihiCh is usually the best of 
1-dimensional minimization methods, to find the appropriate 
step length a*. The interior penalty function has been used 
to transform the problem into an eqxiivalent unconstrained 
optimization problem. 

In interior poialty function method, a 0 function is 
constructed by augmenting a penalty term to the objective 
function F(X) as given by Eq, (1.3) 
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By tlie repeated vmconstrained minimization of the 
J2f-fvinetion for a sequence of decreasing values of the penalty 
parameter (k* 1, 2 , 3 ...)# the solution converges to that 
of the original problem stated in Eq. (1.2). Hence ^ penalty 
function methods are also known as sequential unconstrained 
minimization techniques (SUMT). The penalty term is chosen in 
such a way that its value will be small at points away from 
the constraint boundaries and will tend to infinity as the 
constraint boundaries are approached. Thus, once the uncons- 
trained minimization of |2f-f unction is started with any feasible 
solution, the subsequent solutions generated will always lie 
within the feasible domain, since the constraint boundaries 
act as barriers during the minimization process. The reason 
for not carrying out the minimization of the {2t-function in 
only one stage using a small value of penalty parameter r, 
so that the minimum of 0(X, r) tends to the minimum of F(X), 
is that 0 will be a very complex function away from the optimum 
and the algorithm will not get initiated if r is small. 

However, near the optimum point, the (2f-f unction is well 
behaved and permits computaticais to be carried out with a 
small value of r. Therefore, minimization is carried out 
as a sequence of unconstrained minimizations with decreasing 
values of r and the final design variable vector of the 
previous cycle is taken as the initial design variable vector 
for the next cycle.. 



18 


1.5 Need for Optimal Design Tables 

Optimal designs of structures with limited parametric 
studies result in not only economical and efficient designs 
but also in better understanding of the response of the 
structure and in identifying the factors which largely influ- 
ence the cost or the efficiency of the structure. Ihus^ the 
incorporation of optimization concepts results in both economy 
and improved techniques. 

In the field of structxiral design, optimization has 
essentially remained a pursuit for research (Goble and Moses, 
1975). No doubt, econcxny and efficiency can only be ensured 
through optimization techniques. But, an average structural 
engineer, particularly in India, is neither farailier with 
optimization methods nor can afford to pay for the necessary 
computer software and time to carry out optimization studies. 

Therefore, for the structures considered in this 
stxady, parametric studies have been made and results are 
presented in the tabular form. In order to enable the design 
engineer to directly pick up the optimal design, optimization 
studies have been carried out within the framework of codes 
of practice. Xs far as possible, the specifications of 
Indian standard codes have been followed for carrying out 
those designs. However, in certain cases of discrepancies 
and recent developments which are not yet incorporated in 
Indian codes, BS 5337 and CP 110 are followed. 



1.6 Organisation of the Thesis 


The thesis is divided into seven chapters of which 
the first one gives the brief introduction of various analysis 
and design philosophies, optimization techniques besides iden- 
tifying the scope of the present study. 

In Chapter 2, a selective review of literature related 
to scope of the thesis is made. Based on literature review a 
need for further investigation on certain aspects has been 
identified. 

Optimal limit state designs of Type I, Type II, and 
Type III water towers are dealt with in Chapters 3, 4 and 5 
respectively. Chapter 3 constitutes the main body of the 
thesis which deals in detail the force method of elastic 
analysis for mxiltishell stiructures limit analysis of the coni- 
cal reservoir, and all other design aspects. The general 
layout of each of these chapters has been, a brief discussion 
of the design considerations followed by the methods of 
analysis, formulaticxi of the optimization problan, presenta- 
tion of optimal designs in the tabular form, discussion of 
the results of parametric study and an eacample to illustrate 
the method of using the design tables. 

Chapter 6 deals with computational aspects in which 
the details of the optimization technique used in this study 
are discussed. Certain important observations made during 
this study, related to computation, are hi^lighted. 

Finally, summary and the ccwiclusions resulting from 
the study along with suggestions for future work are presented 
in Chapter 7. 



CHAPTER 2 


LITERATURE REVIEW 

Hie present study enccxnpasses many fields such as the 
response of reinforced concrete structures to service loads / 
limit analysis of structures/ limit state philosophy and opti- 
mization. Literature review on some of these fields has been 
well documented and available even in some textboo]cs (Save 
and Massonnet/ 1972? Green and Perkins / 1980). Hence only a 
selective review of literature/ on the elastic analysis of 
shells, limit state philosophy, and optimization technique and 
its application to reinforced concrete water towers will be 
made in this chapter. A brief review of literature for limit 
analysis of reinforced concrete shells and cracking phenomenon 
of reinforced concrete under service loads has been given at 
appropriate places. 

2.1 Elastic Analysis of Shells 

The shell theories based on linear elasticity concept 
are most commonly used. These theories adequately predict 
stresses and deformations for shells ejdiibiting small defor- 
mations, The nonlinear theory of elasticity forms the basis 
for finite-deflecticxi and stability theories of shells, 
Large-deflectica:i theories are oftai required when dealing 
with shallow shells, highly elastic membranes and bucking 
problems. The nonlinear ^ell equations are considerably 
more difficult to solve and for this reason are more limited 
in use. 



Love (1944) was the first investigator who presented 
fir St -order-approximation shell theory based on classical 
elasticity and gave simplified strain-displacement relations 
and consequently the constitutive relations. 

Second-order-approximation shell theories retain the 
t/R terms, where ‘t* is the shell thickness and ' R' is the 
least radius of curvature of the middle surface, in comparison 
with unity in the stress resultant equations and in the strain 
displacement relations. Several authors (Flugge, 1960; 

Kemjner, 1955; Novozhilov, 1959 and Kraus, 1967) have contri- 
buted to the development of these theories. 

Foregoing theories take into account the flexural 
behaviour of shells and generally referred to as 'bending* 
theories of shell . Using these theories solutions for special 
cases have been found by many investigators (Pfluger, 1961; 
Salvador!, 1955; Baltrukonis, 1959; Galletly, 1955 and 1960 
and Baker and Cline, 1962). 

A special case of bending theory of shells in which, 
in the study of equilibrium of a shell, all moment expressicns 
are neglected, is known as the man(±)rane theory of shells. 

First contribution to membrane theory were furnished by Ijame 
and Clapeyron early in nineteenth centiury. They considered 
symmetric loading on shells of revolution. A number of authors 
(Novozhilov, 1959; Goldenveiser , 1961 and Timoshenko and 
Krieger, 1959) have contributed useful material in the deve- 
lopment of membrane theory of shells. Using this theory, 
soluticms to innumerable special cases of practical importance 



have been worked out by several investigators (Pfluger^ 1961; 
Flugge^ 1960; Salvador!, 1955 and Haas, 1962). 

Timoshoiko (1940) gave analytical membrane solutions 
for axisymmetrically loaded orthotropic shells of revolution. 
Baker (1964) derived governing differential equation for 
ortho tropic composite cylinders subjected to axisynroetric 
loads and obtained results for some special cases of practical 
importance. 

Gerard (1962) presented solutions for orthotropic 
cylinders under axial compression using linear stability 
theory. Becker and Gerard (1962) derived small deflection 
theory for general instability of orthotropic circular cylin- 
drical shells for external pressure, torsion and axial compre- 
ssion. Experimental data were compared with the theoretical 
results and was found to agree reasonably well. 

Elastic analysis and design of reinforced concrete 
water tanlcs, with relevant literature, has been given in some 
textbooks (Green and Perkin, 1980; Gray and Manning, 1973; 
Krishna and Jain, 1980; Baikov, 1978), 

Development of more exact theoretical equations does 
not necessarily assist in the solutiOTi of practical shell 
problems, because the resulting theoretical equations can 
only be solved with great difficulty, and even then only for 
special cases. Ihe experimoital approach is also limited 
because data are not available for every special case. 
Practical difficulties in both theory and e3q>eriinent have 
led to the developm^t and application of applied engineering 



methods for the analysis of shells, While these methods are 
approximate and are valid only under specific conditions, 
they generally are very useful and give good accuracy for the 
analysis of practical engineering shell structures. 

?irya (1969-70) obtained, using bending theory of 
shells, numerical data and presented in a convenient tabular 
form so that the continuity effect at the junctions of a multi- 
shell structure could be carried out with ease and accuracy. 

The use of these data was explained through the analysis of 
an Intze tank. Kameswara Rao (1982) also presented nxjunerical 
coefficients for rapid calculations of continxiity effect in 
Intze tank. 

Baker et al. (1972) have provided a very good collec- 
tion of membrane and bending solutions which are usefvil for 
the analysis of axisymmetrically loaded shells of revolution 
Qicountered in engineering practice. The force method of 
analysis specially suited for such shell structures, also has 
been dealt in detail. 

Kalwar et al. (1983) described the force method of 
analysis for axisyrametric multishell structures and applied 
it to the analysis of a conical water tank supported cai 
cylindrical shaft. J^didam et al. (1985) have also used this 
method for the analysis. 

2.2 Limit State Philosophy 

ThouiQptit was given to limit state philosophy in Soviet 
Union as early as 1930 and the limit state approach was embo- 
died in the Russian Code in 1954 (Zalesov, 1973). 
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The question of what constitutes an acceptable risk 
or probability of failure is very controversial. Society has 
long demanded and still doeS/ a degree of safety in structures 
which is out of proportion to the risk it accepts in other 
areas of life. As this concept of safety is deep rooted in 
structural design, the notion of probability of failvxre, even 
if it is of the order of 10 , is repulsive to the majority of 

present day engineers. This is reflected by Baker's statement 
(1966) that except for absolutely unforeseeable events, a zero 
probability of failure has to be demanded from the designer. 
But, to quote Prof. Freud enthal, 'The difference between safe 
and unsafe design is in the degree of risk considered accep- 
table, not in the delusion that such a risk can be completely 
eliminated' , 

For the first time a probabilistic approach to the 
structural safety, taking into account expected variations in 
load and strength, was introduced in the structural design in 
the aircraft industry in the late 1940s. 

Following the attempts by Freudenthal (1947, 1948) 
to derive the safety factor from statistical variations of 
the design parameters, concept of safety was extensively 
considered by Pugsley (1951) and Johnson (1953). 

A composite international body of research workers 
and practising designers known as Exoropean Committee of 
Concrete, CEB, was set up in 1953 with the object of providing 
a sound idiilosophy of design, together with recommendations 
on the detailed aspects of design based upon ejqserimental 
and theoretical research. With the inception of CEB much 



more serious attention began to be devoted to the entire 
subject of structvuTal design in reinforced concrete. Ihe 
limit state philosophy has been adopted and elaborated by the 
CEB and formed the basis of that Committee' s recommendation 
(CEB, 1964). 

Rowe (1965) forwarded the concept of limit state 
design and stated that this philosophy is simply a restatement, 
in a logical and rational manner, of the essential concepts 
of design which has been obtained for many years. 

Cornell (1969) proposed a code format for necessary 
safety factor in the struct\iral design. He suggested to 
incorporate the measures of uncertainty in predicting the 
strength and applied force and a measure of degree of safety 
in determining the safety factor, 

Kemp (1973) gave a number of valid reasons why class- 
ical reliability theory cannot be applied, in its purest form, 
to the design of reinforced concrete structures. 

According to the International Standards Organization 
(TC98, 1973), it is not generally required that structures 
should be designed to withstand effects of exceptional events 
like wars. In case of certain other events, like accidents 
or earthquakes, when their frequency is ill-defined, the 
designer should ensure that the risks associated with such 
events are limited. 

Ang and Cornell (1974) also pointed out that proba- 
bility of an adverse event or ' failure* is virtually unavoi- 
dable and proposed that the failure should be interpreted 



with respect to some predefined limit state. Depending upon 
the limit state under consideration/ the concept of fallxire 
probability can be applied to both the safety and the perfor- 
mance of structure. Failure probability need not/ however/ 
appear explicitly in the routine design equations. 

Ellingwood and Ang (1974) developed the methods to 
illustrate the quantitative analysis of uncertainties and 
showed the effect of these uncertainties on the level of risk. 
At the same time risk associated with existing design proce- 
dures with reference to reinforced concrete were evaluated 
which formed the basis for consistait code development. 

Following the attempt by Ravindra et al. (1974) to 
illustrate reliability-based designs of simple structural 
elajnents, Moses (1974) considered the reliability of struc- 
txiral system to safety. 

Attempts have been made by Galambos et al. (1982) 
and Ellingwood et al, (1982) to describe the methodologies to 
develop probability-based load factors and loading combina- 
tions for use with common construction material and technolo- 
gies and to assess structural resistance. 

Ellingwood (1982) examined the relative advantages 
of a number of alternate safety checking formats with regard 
to their ability to provide uniform reliability in limit 
state design, 

Ellingwood et al. (1983) suggested that explicit 
provisions to mitigate the effects of atariormal loads should 
be a part of building codes and standards. 



since code formulation is an evolutionary process, 
the new provisions seldom reflect abrupt changes from estab- 
lished ones.. The International Standards Organisation, 
aiming at unification of different methods of structural 
calculations and ensiiring safety of structures, has recotmi- 
ended a semi-probabilistic limit state method. On the basis 
of these recommendations, codes of practice for design of 
reinforced concrete structures are being revised to conform 
to the limit state philosophy, BS 5337, CP 110, IS; 456 are 
amcsng those in which the limit state theory is incorporated 
and are referred to in the present work. 

As the present codes of practice are semi-probabi- 
listic in natvLce, the futxare codes of practice have to be 
modified drastically to incorporate the data that are likely 
to be available in the near future, 

2.3 Optimization Technique and Reinforced Concrete Water 
Towers 

2.3.1 General 

The object of a structural design is to achieve not 
only a safe and a serviceable structure but also the • best* 
and 'economical* one. The methods of achieving optimum 
designs are primarily concerned with arriving at the 'best* 
structure according to some criterion of structural efficiency. 
Usual criteria considered for such optimization are the 
minimum wei^t/cost, or the maximum stiffness etc. 

Excellent survey made by Wastiutynski and Brandt 
(1963) and Sheu and Prager (1968) provide a complete 


historical development of the methods of structural optimi- 
zation , 


Fletcher and Powell (1963) gave a very stable rapidly 
converging method for unconstrained nonlinear programming 
problems of optimization. Whereas Piacco and McCormick (1968) 
developed sequential unconstrained minimization tecdanique for 
the solution of constrained nonlinear minimization problem by 
transforming it into an equivalent unconstrained problem of 
minimization , 

Fox (1971) and Rao (1978) have dealt with the various 
methods of optimization and givon relevant literature survey 
in the field, 

2,3,2 Application to reinforced concrete water towers 

In this section, literature mainly related to appli- 
cation of optimization techniques to reinforced concrete water 
towers has been covered, 

Meldiers and Rozvany (1970) solved the problem of 
minimizing the reinforcement in cylindrical reinforced concrete 
water tanks of a given thickness using Prager-Shield opti- 
mality conditicxi. The most important aspect of cracking was 
ignored by them in the design, 

Nielsen (1974) gave a procedure to find the required 
area of reinforconent at any point of a shell surface using 
the plastic design method. Thickness of the ^ell and the 
distribution of the internal forces were assumed to be known. 
Depending on the relative magnitxide of the bending moments. 



twisting moments, and the membrane forces, the optimal pattern 
of the reinforc«T>ent was chosen, 

Wilby (1977a) in his extensive literature sxirvey on 
str\ictural analysis of reinforced concrete tanks has concluded 
that very little work has been done on the optimization of 
tanks using either elastic or plastic models. 

Selvanathan (1978) has considered the use of different 
number of columns and spacing of bracings to arrive at an 
optimal design of the Intze tank for two capacities using 
nonlinear programming. The total cost of stsructure, which 
includes cost of tank, staging and foundation, has been chosai 
as the objective function, The working stress method is 
adopted for the design. 

Bond ( 1979 ) proposed a method , which combines iso- 
parametric finite element analysis with a nonlinear optimi- 
zation technique, to search for structural forms resulting in 
minimum cost reinforced concrete structures. One of the 
examples considered is that of a conical water tank. The top 
diameter and thicknesses at top and bottom have been chosen 
as design variables. The objective function includes the 
cost of tank, support and the top dome. Alternative method 
of design, which is given in BS 5337, has been used to govern 
the resistance to cracking. Two Important parameters, one 
being the reinforcement and its distribution and the other 
being the diameter of the supporting structure, are not taken 
as variables in the optimal design process. 

Subramanyam (1981) and Adidam and Subramanyam (1982) 
investigated the optimal design of surface water tanks using 
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the limit state design philosojhy, Radius of the tank, 
thickness of the tank wall and the reinforcement areas in the 
two directions are chosen as design variables. Cost of 
finished concrete, reinforcement and formwork have been taken 
as the objective function. The permissible surface crack 
width has been taken afe equal to 0,2 mm. which seems to be 
higher for water retaining structxures and the affect of 
seismic forces is ignored in the design. Resistance to 
cracking due to hoop tension is ascertained by limiting the 
stress in reinforc«t\ent only without assigning any upper limit 
for direct tensile stress in concrete which may result in 
excessive tension in concrete, thereby, causing wide cracks 
in some cases. 

Kalwar et al* (1984) considered the optimal design 
of reinforced concrete water tower consisting of a conical 
tank supported on cylindrical shaft using limit state design 
philosophy «ind a nonlinear optimization technique. The 
diameter of the supporting ^aft, ratio of the top diameter 
to bottom diameter of the conical reservoir, thicknesses of 
various elen^nts and the reinforcements in these elements 
with practicable logical distribution of reinforcement are 
taken as design variables. The objective function includes 
the costs of materials, formwork and labour for superstruc-s 
t\ire. The procedure given in IS; 3370, with a little modi- 
fication with regard to permissible stress in steel, has 
been used to provide resistance to cracking. The effect of 
escalation in unit costs, of materials and formwork, on 



opitimal configuration and the design variables is also inves- 
tigated. 

The following studies do not make use of any optimi- 
zation techniques in the strict sense , but consider ways of 
minimizing cost or weight of reinforced concrete water tanks, 

Wilby (1977b, 1978) has considered the cost minimiza- 
tion of polygonal tanks, in stages. In the first stage, for 
a given capacity of the tank, the surface area is minimized 
with a view to minimize the cost of the shuttering. Then the 
thickness of the rectangular plates and reinforcement in them 
are determined from individxial yield-line collapse mechanisms.’ 

Jain et al, (1979) have carried out extensive para- 
metric studies to arrive at optimal designs of Intze tanks. 
Various capacities, staging heights, bearing capacities of 
soils and lateral forces due to wind or earthquake are consi- 
dered for the analysis. Based on these studies, a number of 
equati<»is for rapid estimation of cost and materials have been 
presented. For computing the seismic forces the seismic 
coefficient method has been used instead of the response 
spectrum method. Nothing has been said about the final 
design variables for any of the cases considered and about 
the procedure adopted to arrive at the optimal designs. 

K^meswara Rao and Raghavan (1982) considered the 
optimal design of Intze tanks on shaft. In this study attempt 
has been made to optimize the weight of the reservoir portion 
only. This has resulted in a large diameter of the supporting 
shaft giving an unpleasant cc«ifiguration of the overall struc- 
ture which may not be economical in the true sense of the 
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word. Moreover# wind forces and the continuity effect at 
joints have been ignored in the design. 

2.4 Concluding Remark 

It is obvious# from the foregoing literature review# 
that application of limit analysis to the design of water 
tanks has hardly been attempted. Optimal design of water 
towers# with the exception of the Intze tank# has received 
very little attention although considerable amount of money 
is spent on the construction of water towers. In a developing 
country like India where a large number of water towers will 
have to be constructed to meet the requirement# optimal design 
of such structures will be very meaningful. 

In view of the foregoing considerations optimal 
design of reinforced concrete water towers has been considered 
in the current study. The emphasis has been to obtain prac- 
tically usable optimal cost designs# conforming to the rele- 
vant codes of practice# which will remain safe and serviceable 
throughout their design life- 



CHAPTER 3 


OPTIMAL DESIGN OF TYPE I WATER TOWERS 


3.1 General 

Amc«ig the various types of reinforced concrete water 
reservoirs or tanks, cylindrical, conical and Intze tanks are 
the most ccanmonly used. Reservoirs with the conical walls 
have a rainimtiin. cross-sectional area at the joints with the 
supporting towers, which makes it possible to efficiently 
establish the dimension of the supporting towers. Moreover, 
conical tanks supported on cylindrical shafts are aesthetic- 
ally pleasing ones and are preferred because of their better 
architectural expressivoiess. In. this chapter optimal design 
of Type I water towers (Pigxire 1,1) is presented. 

3.2 Design Loads 

Various loads to be considered for the design of a 
water tower are dead load, hydrostatic pressure loading, 
weight of water, wind, and seismic loads. 

The values of wind and seismic loads depend upon the 
geographical locaticm of the place where the tower is going 
to be situated. India is divided into three wind and five 
seismic zones depending upon the intensity of wind pressure 
and seisttnicity. 

3.2.1 Wind loads 

IS; 875 specifies the values of basic wind pressures 
for the various wind zones in India. These basic wind 



pressures are nothing but the equivalent static wind press- 
ures which take into account the dynamic head pulsaticai due 
to wind gusts. Ihe basic wind pressures for these zones are 
given as -under 


Wind zone 
Zone I 
Zone II 
Zone III 


Basic wind pressure * 
1^0 kN/m^ 

1.5 

2.0 WSI/ra^ 


♦■For heights upto 30 m above the main retarding surface. For 
heights more than 30 m IS: 875 is to be referred. 


IS: 87.5 also specifies^ for structures of various 
plan shapes other than rectangular ones^ the external wind 
pressure acting on the projected area in the plane perpendi- 
cular to the wind, as the product of basic wind pressure and 
the shape factor. For structures having circular plan shapes, 
the specified value of shape factor is 0.7. Hence, in calcu- 
lating wind loads the Shape factor is taken as 0,7.. 


3,2.2 Seismic loads 

For the pxirpose of determining the seismic forces, 
the country is classified into five zones, IS: 1893 specifies 

■the values of basic seismic coefficients (a ) and seismic 

o 

zone factors (P^) for computing the design values of hori- 
zontal seismic coefficiaat Depending on the problem., 

one of the following two roe-thods may be used for computing 
the seismic forces 


(a) Seismic coefficient method, and 

(b) Response spectrum method. 

Tile Table 3,1 gives the values of and for different 
zones , 


Table 3,1 Values of basic seismic coefficients and seismic 
zone factors for different zones 


"I 

Zone No, 

— • 

[ Method 


Seismic coefficient method 

1 

1 

1 

.^. 4 — 

Respcxise spectrxam method 

. - 

Basic horizontal seismic 

coefficient, a„ 
o 

1 

1 

1 

1 

t 

Seismic zone factor, 

F 

O 

V 

0.08 


0,40 

IV 

0.05 


0.25 

III 

0.04 


0.20 

II 

0.02 


0.10 

I 

0.01 


0.05 


The values of horizontal seismic coefficients (a, ), 

n 

in the seismic coefficient and response spectrum methods, are 
given by the following ej^essions: 

(a) In seismic coefficient method 

\ » pi (3,1) 

(b) In response spectrum method 

S 

» P I F -S. 
h o g 


(3.2) 




where 


p « a coefficient depending upon the soil -fotandat ion 
system, 

I *= a coefficient depending upon the importance of 
the structure, 

a© = basic seismic coefficient as given in Table 3.1, 
P© seismic zone factor for average acceleration 
spectra as given in Table 3.1, and 
f = average acceleration coefficient as read from 
Figure 3.1 for appropriate natural period and 
damping of the structure. 

In the present study, for finding the wind and 
seismic forces, the procedures specified in Sections 3.Z-1 
and 3.2,2 have been adopted. For finding the value of average 
acceleration coeftigient, the damping of reinforced concrete 
structure is taben as 5%. 

3.2.3 Loading combinations and partial safety factors for 
loads 

It is not always necessary to analyse a structure for 
all the possible loading combinations for arriving at the 
governing design forces in different parts of a structure, 
Therefore# it is usefxil to look critically at the various 
possible loading combinations, to decide which of these 
combinations are most likely to govern the design of a 
particular portlcan of the structure. Based upon the preli- 
minary studies conducted by the author and the experience 
of the several designers and investigators in India , Arya 





(1969^ 1970), Jain and Singh (1977), Jain et al . (1979), 
Kameswara Rao (1982) and Karneswara Rao and Ra gha van ( 1982) it 
is recognised that, for elevated water tanks under Indian 
environmental conditions, the governing loading combination, 
for the reservoir portion is the ' dead load + load due to 
hydrostatic pressure' under tank full ccxidition. Whereas, the 
loading combinaticxi of dead load with wind or seismic forces 
under full or a:npty conditions of the reservoir usually 
govern the design of supporting structure. Depending upon 
the loading combination and the limit state under consider- 
ation, IS: 456 specifies partial safety factors for loads,. 

In the li< 3 ht of the foregoing discussions, the loading combi- 
nations and the partial safety factors considered for compu- 
ting the design loads for various portions of the water tower 
are given in Table 3.2. 

3,3 Elastic Analysis 

In order to predict the response of the structure at 
service loads, elastic analysis is carried out. Water towers 
considered in the presaat study can be regarded as m\iltishell 
structural systems because they are composed of cylindrical , 
conical and spherical shells and circular rings. These towers 
are axisymmetric in nature, A shell , in general , is a stati- 
cally indeterminate structure. Hitherto, thin ^ell theories 
based on linear elasticity concepts are most ccxtimonly iised. 

The exact solution for a shell structure is very difficult 
to obtain even for an axisymmetric case. As a result, design 
engineers derived solutions based upon simple membrane theory. 
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Table 3.2 Load combinations considered and partial safety 
factor used for loads 


1 

Portion of 

‘^*1 

Load combination 

"■ ■' ■ - ■■ ■ ~ " ' ' ' " ^ 

Partial safety factor 

the tower 
under consi- 
deration 

considered 

- - - 

Limit states 
of service- 
ability 

1 

; Limit state 

{ of collapse 

1 

« 

-—1 



DL LL WL/EL DL LL WL/EL 

Reservoir 

portion 

Dead load (DL) + 
water pressure under 
tank full condition 
(LL) 

1.0 1.0 

1.5 1.5 

Supporting 

shaft 

(a) Dead load + 

weight of water 
lander tank full 
condition (LL) 

1.0 1.0 

1.5 1.5 


(b) Dead load + 
wind/seismic 
forces \inder 
onpty tank 
condition 
(WL/EL) 

1.0 - 1.0 

1.5 - 1.5 

or 

0.9* 


(c) Dead load + 
wind/seismic 
forces under 
tank full 
condition 

1.0 0.8 0.8 

1.2 1,2 1.2 


*This value is considered when stability against overturning 
or stress reversal is critical. 
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However / the solutions so obtained are not found to be satis- 
factory for most of the practical engineering problems. On 
the other hand, the bending theosry of shells is more general 
and exact but vin fortunately it is hi^ly ccxnplex. However, 
for the class of rotationally syiwnetric shells subjected to 
axisytranetric loading, both the theories can be combined to 
arrive at a near-exact solution. The 'force method' is one 
such powerful technique. This method of analysis is very 
well described by Baker et al. (1972) . Kalwar et al* (1983) 
have ejqplained this technique with modified sign conventions, 
better suited for computer aided analysis, and the same is 
employed here for finding the design forces in the reservoir 
portion of the tower. The analysis of the supporting struc- 
ture is carried out by considering the ^aft as a cantilever 
fixed at the base, subjected to wind or seismic forces \ander 
full and empty conditions of the reservoir. 

3.3.1 Elastic analysis for finding the design forces in 
the reservoir portiam 

As discussed in Section 3,2.3, the governing loading 
combination for the reservoir portion of the water tower is 
the dead load + load due to hydrostatic pressure under tank 
full conditicxi, which conforms to the analysis of a rotati- 
onally symmetric multishell structure subjected to axisy- 
mmetric loading. Hence, the * force method* of shell analysis, 
briefly describoi in the next section, is used to analyse 
the tower for finding the design forces in the reservoir 
portion. 



It is worth bringing into light a few facts prior to 
the introduction of the force method of analysis of axisyrame- 
tric shells. When the two solutions, one obtained by simple 
membrane theory and the other by exact berxiing theory, are 
compared, a few useful observations can be made. 

(a) The stresses and deformations obtained by the two methods 
are almost idaitical for all locations of the shell, with 
the exception of a narrow strip cxi the shell surface 
which is adjacent to the boundary or a junction. This 
narrow strip is generally no wider than ''^iRt, where R is 
the radius and t is the thickness of the shell in case 
of spherical and cylindrical shells but for a conical 
shell R is the radius of the circle at the larger end. 

(bj Except for the strip along the boundary or a juncticxj, 
all balding mcanents, twisting moments and transverse 
shears are negligible, this causes the entire solution 
to be practically identical to that of the membrane 
solution. 

(c) In case of the built-in edges, the disturbances along 
the edge are significant. However, the local bending 
and shear decrease rapidly along the meridian and may 
become negligible outside the narrow strip described 
in (a). 

3,3. 1,1 Force method 

The force method of solution for axisymmetric cases 
CTn be explained with the aid of a simple example of a spherical 



shell with built-in edges subjected to an axisyrometric loading 
as shown in Figtare 3, 2(a). 

Since the baiding and membrane theories give practi- 
cally the same result# except for a strip adjacent to the 
boundary# the simple membrane theory can be applied to obtain 
the solutions for the cases shown in Figure 3.2(b). The 
Figure 3, 2(c) shows the corrective forces# moments and hori- 
zontal shears respectively# at the edge which will bring the 
displaced edge of the shell into the position prescribed by 
boundary conditions. Ihe amount of edge corrective forces# 

M and H# will depend upon the magnitude of the deformations 
caused by the primary loading. The final solution can be 
obtained Iny superposing the membrane solutions for primary 
loadings and the solutions for corrective loadings obtained 
by bending theory. This method of solving the problem is also 
loiown as the unit-edge- force method. This method is quite 
similar to the force method used to analyse the discrete 
structural systems. 

The force method explained for a simple element# can 
be extended to multishell problems in a similar manner by 
blowing the multishell structure into simple shell elements# 
for vdiich both primary and secondary solutions are readily 
available# and then impose the equilibrium and compatibility 
conditicxis at the junctions to arrive at the solution. 

The Figure 3.3 shows the water tower mder consider- 
ation# whereas tJie Figure 3,4 shows the same# exploded into 
simple shell elemoits, with corrective loadings at various 
Junctions shown along their positive direction. It should be 
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noted that at junction 3# the elements are not connected to 
each other and hence no corrective force appears at this 
junction. However, even if the elements meeting at this 
junction are connected rigidly, the amount of secondaiTi^ forces 
obtained by actual analysis is found to be negligible. This 
is because of the fact that the deformations of these elements 
at this junction, due to primary loading is very small. More- 
over, due to the large size of all the elements, the deforma- 
tions at one end due to secondary forces acting at the other 
corresponding end of these elem«its is zero. Therefore, 
treating the shell elements meeting at this junction as 
connected or otherwise will not cause any error in the solu- 
tion of the overall problem. This simplification reduces the 
size of the problem and also helps considerably in saving the 
computational time, pairticularly in the framework of optimum 
design , 

The sign convention for corrective edge-loadings and 
deformations at junctions, used for solution by force method, 
can be stated as under. 

Right-hand side: 

All moments in the clockwise direction, 
horizontal loads directed outwardly, 
outward displacements, and 

clockwise rotations at all the junctions, are considered to 
be positive. 
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Left-hand side: 

All moments in the anti-clockwise direction, 
horizontal loads directed outwardly , 
outward displacements, and 

anti-clockwise rotaticais, at all the juncticMis, are considered 
positive, 

3.3, 1.2 Compatibility of deformations and equilibrium 
conditions 

nie compatibility of deformations and equilibrixim 
conditions to be satisfied at various junctions are as follows; 

Junction 1; 

At this junction three elements, top dome, ring beam 
and conical shell, are rigidly connected to eadh other, Ihe 
compatibility of deformations leads us to the fact that rota- 
tion and displacement of each element, meeting at this junc- 
tion, due to the combined effect of various primary and secon- 
dary loadings ^ould be the same. Hence the following equa- 
tions can be written* 


^CNl ^ gTDl 

•CNl ^ .CNl .TDl 


(3.3) 


or 



and 

gCNl gCNl ^^Bl 

P S P 

+ 6g 

(3.4) 


vidiere 6^^^, 6*^^ and 6®^ are the horizontal displac«iients of 
the cOTiical shell, the top dome and the ring beam respectively 
at the junction 1; whereas 6^^/ and and 6^^, 6^^ 

and are the horizontal displacements of these elements 
due to primary and secondary loadings respectively at the 
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The Eqs , ( 3 ..3 ) 

and 

(3. 

4) can be rewritten as 


1 - 

5 TOI 


gTOl 

- 6 ™^ 

(3.5) 

s 

s 


P 

P 

and 6 CN 1 . 

5 B 1 


5 BI 


(3.6) 

s 

s 


P 

P 

similar compatibility 

equations for 

rotation at junction 

1 can 

be written as 






nCNl ^ 
S 



oTDl 

Pp 

«CN1 
- Pp 

(3.7) 

Pfi - 


sr 

ftBl 

Kp 

. 6^1 

Kp 

(3.8) 


where rota- 

tions of the conical shell, the top dome and the ring beam at 
the junction 1 due to primary and secondary loadings respec- 
tively . 

It is important to note that, in this particular case, 
the size of each element is such that the disturbances caused 
by the corrective loadings at one end will die at a short 
distance and will not produce any disturbance at the other «id. 
Hence, the secondary defoirmations at one edge of an elment 
will be the function of corrective forces acting at that end 
only and can be expressed as 


fiTDl 

S 

,TD1 
“ *11 

XI 

^ ,TD1 
^12 

X2 

(3 

.9) 

qTdi 

.TDl 

f 

21 

XI 

4. 

^22 

X2 

(3 

. 10 ) 


The iqs. (3.9) and (3.10) can be written in a convenient form 
using matrix notations as \jnder 
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where 


and 


atid the 


and 


,TD1 




,TD1 




— 


w- — 


.TDl 

*11 

.TDl 

*12 


XI 


.TDl 
_ 21 

.TDl 
22 _ 


X2 


(3,11) 


-TDl 


11 


,TD1 

■12 


: horizontal displacement of the top dome at 
the jimction 1 due to a positive unit hori- 
zontal force per unit length acting at that 
junction , 

horizontal displacement of the top dome at 
the junction 1 due to a positive unit moment 
per \init length acting at that junction, 
rotaticn of the top dome at the junction 1 
due to a positive \init horizontal force per 
unit length acting at that junction, 
rotation of the top dome at the junction 1 
due to a positive unit moment per unit length 
acting at that junction. 

Similar expressions can be written for the ring beam 
conical shell at the junction 1 as under 


pTDl 

■21 


.TDl 

'22 


6 

^ i 
1 

p! 


B1 


.B1 -Bl 

S 

tx 

11 12 

B1 


.Bl .Bl 

S 


21 *22 



r “■ 

CNl ' 


CNl CNl 

S ‘ 


11 12 


s 


CNl 


CNl CNl 

S j 


_ 21 22 _ 


X3 

X4 

X5 

X6 


(3.12) 


(3.13) 


Ihe equilibrium conditions will yield two more equa- 


tions for the juncticm 1 as follows 
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XI + X3 + X5 = 0 (3.14) 

X2 + X4 + X6 = 0 (3.15) 

The relaticxis (3.11) through (3.13)^ the compatibi- 
lity conditions (3.5) through (3,8) together with the equili- 
brium Eqs. (3.14) and (3.15) can be arranged as a set of six 
linear algebraic equaticais in six unKnowns (Xl to X6) as i;under 


.TDl 

'^11 

.TDl 

"*12 

0 

0 

.CNl 

^11 

XNl 


XI 


■gTDl 

P 

- 


.TDl 

21 

,TD1 

22 

0 

0 

XNl 

^21 

.CNl 

22 


X2 



- 


0 

0 

^11 

31 

“^12 

.CNl 

11 

fCNl 

^12 


X3 


gBl 

- 

, CN 1 

0 

0 

rH r-i 

m 

1 

-f®^ 

*^22 

.CNl 

^21 

fCNl 

^22 


X4 


oBl 

Pp 

- 


1 

0 

1 

0 

1 

0 


X5 



0 


0 

1 

0 

1 

0 

1 


X6 



0 









- 






(3.16) 


Junction 2; 

At this junction four elemaits, the ring beam^ the 
conical shelly the supporting cylindrical shell and the inner 
cylindrical shell are rigidly connected to each other. Procee- 
ding in the same way, the compatibility and equilibrium condi- 
tions can be expressed by a relation similar to that givon 
by Eq. (3.16) for the junction 2. Ihis relation will represent 
another set of eight linear equations with eight unknowns 
(X7 to X14) as follows. 
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cCN2 


where SZ‘^% 5®^, 6;" and 8;“ and Pp’‘, and p; 


S2 

'P 


CY2 
•p 


,CN2 


,B2 „S2 


CY2 


are the displacements and rotations at junction 2, of conical 
shell, ring beam, supporting shaft and inner cylindrical 
shell resi>ectivelY , due to primary loadings 
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.CN2 

■21 


.CN2 

'22 


= horizontal displacement of the conical shell at the 
jvmction 2 due to a positive unit horizontal force 
per \init length acting at that jtinction, 

= horizontal displacement of the conical shell at the 
jmction 2 due to a positive \jnit moment per unit 
length acting at that jxanction, 

= rotation of the conical shell at the jvinction 2 due 
to a positive unit horizontal force per vinit length 
acting at that junction, and 
= rotatictti of the conical shell at the junction 2 due 
to a positive unit moment per unit length acting at 
that junctic»i . 
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Similarly , 

0Y2 '■'CY''? 

£21 ^22 ^*P^®sent the displacements and rotations at 

junction 2 due to unit edge loadings for ring beam, cylin- 
drical support, and the inner cylinder respectively. 


As already discussed earlier, the size of all the 
shell elements in this particular case are such that the 
influence of the corrective loadings at one end do not inter- 
act with those caused by the similar loading at the other 
end. Hence, the two sets of equations, (3,16) and (3,17), 
can be solved independently or combined together to form a 
single set of fourteen linear algel:^|^l^qpati<^j^^j||Eourteen 

« W 


un}<nowns to arrive at the solution. In certain cases, where 
the stresses and deformations at one end are influenced by 
the corrective-edge loadings at the other end, the soluticn 
can be obtained caily by solving the combined set of linear 
relations involving all the unknowns. 

When Eqs, (3.16) and (3.17/) are combined together, 
the res\ilting equation will be of the form as follows 

[f] X = Ap (3.18) 

in which [f ] is the flexibility matrix, 

X is the load vector of the unknown corrective 
edge-loading, and 

Ap is the deformation vector due to primary 
loading. 

It can be seen from Eq. (3.18) that the solution by 
force method requires formulae for obtaining the elements of 
the flexibility matrix and the deformations of the various 
shell elements at the junctions due to primary loadings. 

Most of the expressions for stress resultants and 
deformations at any point of the shell, needed for the present 
work, are documented in the existing literature. A few such 
expressions, for specific cases, are derived in the present 
work. All these are listed in tabular form for various 
primary and edge-loadings. These expressions are given for 
both special and general cases to obtain defoimations at 
junctions only and the stress resultants and deformations 
at any point of the shell respectively. For primary loading 



solutions are obtained by using membrane theory whereas 
bending theory is used for the edge-loading. 


3. 3. 1.3 Special solutions for deformations at jiinctions 


These expressions are nothing but a special case 
of the general load-deformation relations at any point of the 
shell. The Tables 3.4 through 3,10 give the formulee for 
deformations at junctions due to various primary and correc- 
tive edge-loadings. The sign-conventions adopted here is the 
same as mentioned in Section 3. 3. 1.1, In order to put the 
expressions for the bending solution in a convenient form for 
conical and spherical shells^ factors and F^(n) are used. 
These factors are given in the Table 3.3 in which the symbols 
K# L and n for conical and spherical shells are given as 
follows. 
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t as thickness of the shelly 
as Poissonfe ratio, 

E =s Youn^ modulus of elasticity, and 

L, 3 ^, X, L and a are as indicated in Figure 3,5. 


Spherical shell: 
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Table 3.3 (continued) 
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Table 3.4 Open conical shelly special solutions for deformations at junctions due to edge- 
loadings 



For T., D and K see Table 3,3 and section 3. 3. 1.3 
















For Fj and k see Table 3.3 and Section 













Table 3.6 Long cylindrical shell, special solutions for deformations at junctions due to edge 
loadings 



For D and % see Section 3,3.1 































Fig. 3-7 Details of cylindrical 
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Table 3.8 Open conical shell, special solutions for deforma ticms at junctions due to primary 
loadings 
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Table 3.8 (continued) 



note:- • f * for the present case is equal to sin a 


Table 3.9 Open spherical shell# special solutions for deformations at junctions due to primary 
loadings 



T^ble 3,10 Cylindrical shell, special solutions for deformations at junctions due to primary 
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In the presait study both, supporting and the inner 
cylinder, come under the category of long cylinders in which 
case > 5, Therefore, the simplified solutions which are 

meant for long cylinders are given for use. 

Now \ising the solutions given in Tables 3,4 through 
3.7 for H = M = 1, the various elements of the flexibility 
matrix [f] can be found. Similarly Tables 3,8 through 3,10 
can be used to obtain the deformations of shell elements at 
varioiis junctions due to primary loadings and then all the 
elements of deformation vector TJp can be determined. The 
Eq, (3.18) can now be solved to obtain the unknown values of 
corrective edge-loadings (XI to X14). 

After obtaining the valxJies of XI to XI 4, the stresses 
and deformations at any point of the shell can be obtained 
by superposing the membrane solutions for primary loadings 
and the balding solutions for the edge-loadings. For this. 



general formulae are needed for both stresses and deformations 
due to primary and edge-loadings* The next section deals with 
these general solutions. 

3, 3. 1,4 General solutions for stress resultants and 
deformations 

The Tables 3.11 through 3^16 give solutions for 
stress resultants and defomatioos at any point of the shell 
under consideration for various cases of primary and edge- 
loadings. The Figure 3.8 shows an infinitesimal element of an 
axisymmetrically loaded shell of revolution. The element is 
in equilibrium under the indicated loadings. The stress 
resultants and shown in Figure 3,8 are taKen 

as positive and all the general solutions given in Tables 3,11 
through 3,16 conform to this sign convention. The nomencla- 
ture used for different stress resultants and the sign conven- 
tion mentioned in the foregoing can be stated as under; 

= force in the meridional direction^ kW/my positive if 
causing tension, 

Niq « force in the circumferential direction, WN/m, positive 
if causing tension, 

M 0 * moment in the meridional direction, kNra/m, positive if 
causing tension on inside face, 

« moment in the circumferential direction, MSlm/my positive 
if causing tension on inside face, and 
12 ^ = tra^ii 3 ?^ReacigB sHewr ftsree in the meridional direction , 
positive as indicated in Figure 3.8, 




Table 3,11 (continued) 





Fig. 3*13 E dge-loQd i ng s on cy i i nd ncai shell 
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Tabl e 3 . 1 2 ( continued ) 
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Table 3.12 (continued) 



Table 3.12 (continued) 



Table 3.13 Bending solutions for long cylindrical shell - edge-loadings 



[(2 V 2 ) MRKj^ e cos(K. n + 71/4) j/L 


^ble 3 » 1 a { continued ) 



Contd 


Table 3.13 (continued) 



For D and K- see Secticw) 



table 3.14 (continued) 




M«nbrane aolutiom tor cylindrical shell - prlewity loadinas 




The sign convaition for deformations in general solutions 
also remain the same as given in Section 3.3. 1.1. 

For a shell of revolution with axisymmetric loading, 
the twisting moments, membrane inplane shears and the trans- 
verse shear in the circumferential direction are zero. 

3. 3. 1,5 Applicability of membrane and bending solutions to 
reinforced concrete shells 

The reinforced concrete shells, in general, are 
orthotropic in nature. The foregoing discussion is applicable 
to axisyrametrically loaded shells of revolution of homogeneous 
isotropic material. The membrane solutions given for finding 
the value of and can also be used directly for ortho- 
tropic shells because, the membrane is a statically determi- 
nate structure and stress resultants do not depend on the 
stiffness properties. 

For special case of orthotropic shells in which the 
oxtensional and bending stiffnesses in meridional direction 
are equal to the respective stiffnesses in circumferential j 
direction, the expressions described in the previous section 
can be used by merely replacing t and E by their modified 
values as described below s 

[ 

In expressions for deformations obtained by manbrane sol\ition j 
for primary loading, we set 

Et » B(1 - (3.19 

and, in expressicwis fear stress resultants and deformations | 
obtained by bending solutions for edge loading, we set 



t 


and 


E » 


(3.20) 




where B « » Bq »* extensional stiffness of the shell in 

meridional or circumferential direction. 

D * = Dq a= balding stiffness of the shell in meri- 

dicsial or circumferential direction. 

Ihe values of B and D can be found by the following expressions 




®i^i 

777 


and 




^i^i 

2 

1 - 


where i refers to or 6# 

A is the cross-sectional area of 1 m wide strip, and 
I is the moment of inertia of 1 m wide strip. 


It is further important to note here that in rein- 
forced concrete water reservoirs the limit state of craclcing 
governs the design. Therefore/ they are designed such that 
tension in the concrete remains well within its cracking 
limit ensuring entire concrete section to be effective. It 
^ould also be noted that, for the loading combination consi- 
dered for reservoir portion, the supporting shaft remains in 
compression only. Even fear other loading combinations, a 
considerable length of the shaft near junction 2 will remain 
in axial compression throtoghout and the tension produced by 
wind induced ring mooMWit and hoop tension will be well 
within the cra^^cing limit, 

AS the section remains uncracked, the contributiem 
to extensional and bending stiffnesses due to the area of 
reinf orcemarvt is small as compared to that contributed by the! 



area of concrete. Taking this into consideration, the stiff- 
nesses of a shell element in the meridional and circumferen- 
tial directions can be taken equal for all practical purposes 
even if the reinforcement in the two directions is slightly 
different. Moreover, these are not the absolute values of 
stiffnesses of any shell element which affect the result but 
are the relative values of stiffnesses of various shell 
elements interacting at the jxinctions. Ihus, even if the 
contributicxi of reinforcement to stiffnesses is neglected, 
the solution so obtained will be quite accurate. In view of 
the foregoing discussion the stiffnesses in the two directions 
of a shell element turn out to be equal, thus facilitating 
the use of the solutions of homogeneous isotropic shells for 
reinforced concrete reservoirs taking E as the Youngs modulus 
of elasticity of concrete and t as the thickness of the shelly 
Results so obtained are reasonably accurate, in particular 
the hoop tensions at locations which usvially govern the desigi 
of conical tank wall remain unaffected. 

3,3*2 Elastic analysis for finding the design forces in 
the supporting shaft 

Ihe analysis of the supporting structure is carried 
out by considering the shaft as a cantilever fixed at the 
base. The various loading conbinations to be considered for i 
finding the design forces in the shaft are those given in 
Table 3.2. The analysis for the case of dead load + weight j 
of water \inder tahic full condition has been already discussei 
in the previous section, in order to carry out the analy sis 



for the other two load combinations, the design values of 
wind pressures and the seismic forces are required to be found 
out. The design wind pressure is obtained by following the 
procedure given in Section 3,2,1 and using an appropriate 
value of partial safety factor given in Table 3,2, 

For finding the design seismic forces, the response 
spectrum method is used to find the value of horizontal 
seismic coefficient. IS; 1893 specifies a simplified method. 
In this method, elevated tanks are considered to be the 
systems with a single degree of freedom with their mass concen* 
trated at their centres of gravity. The period of free vib- 
ration T, in seconds, is given by 

T « (3.21 

where A =* the static horizontal deflection at the top of 
the tank under a static horizontal force ecjual 
to a weight W acting at centre of gravity of 
tank; and 

g *= acceleratiOTi due to gravity. j 

The val\xe of the weight w, to be used in calculating j 

the deflection A is taken as equal to the dead load of the ! 

tank and one-thiid the weight of the supporting structure 

under empty tank conditicxi. when tank is full, the weight 

of contents is to be added to the weight mder empty conditio 

Using the period T as calculated by Eq, (3,21) and 

damping as 5 percent of the critical value for the concrete 

S : 

Structure, the value of average acceleration coefficient ~ ; 



can be read from average acceleration spectra given in Figure 

3.1. 

In order to find out the value of horizontal seianic 
coefficient using Fq. (3.2), suitable values for p and I 
are to be adopted. For guidance IS: 1893-1975 has specified 
the values for these coefficients. In the present study the 
value of I is taken as 1.5 and the value of P as 1.0 as it 
covers most of the types of foundation under all soil condi- 
tions. 

The lateral force is taken as equal to: 

V 

where * value of horizontal seismic coefficient as 

obtained by Eq, (3.2), and 
W * weight as defined earlier. 

This force is multiplied by appropriate partial 
safety factor to obtain the design lateral force. The design 
lateral force shall be assxmed to act at the centre of gravif 
of the tank horizontally in the plane in which the structure; 
is assumed to oscillate. Ihis is for the purposes of 
carrying out the analysis of the supporting structure. I 

Once the design values of wind and seismic forces 
are obtained, the axial force and the bending moment at any | 
horizontal section of the support can easily be obtained, 

Ihe most stressed section is the one at the top of the j 

foundation. i 

In longitudinal cross -secticns of shell -type 
supports horizcxital loads, particularly wind loads, induce I 



bending moments^ normal/ and transverse forces. In cylind- 
rical supports they may be found as a first approximation by 
analysing elastic rings 1 m high visualized in the support 
between two horizontal cross-sections and equilibrated by the 
horizcwital load acting on the tower at the level where the 
considered ring is located and by the internal tangential 
forces that affect the ring cross-sections. 'Ihe approximate 
value of the induced ring moment, when the wind load is 
distributed over the circumference, is equal to: 

0.33 Wp r^ in KNm/m> 

where Wp = design value of wind pressure at any level in 

kN/m^, and 

r »s mean radius of the shaft in m. 

The hoop force and the shear force come out to be of negli- 
gible value. 

3,4 Design Considerations 
3.4,1 Materials 

The recommended grades of concrete for reinforced 
ccaicrete water tanks, according to IS: 3370 (Part I), are 
M20, M25, etc, but a mix weaker than M20 is not used. In 
view of the early thermal cracking problem associated with 
richer mixes, M20 grade is most ccxnmonly used and is 
preferred in the curramt study also. The minimum quantity 

3 

of cement is kept not less than 3.3 Wl/m in the concrete 
mix. In order to make dense impervious concrete, well graded 



aggregate with adequate water-cement ratio and means of 
thorou^ compaction should be used. It should be recognized 
that common and more serious causes of leaikage, other than 
cracking / are defects such as segregation and honeycombing 
and in particular all joints are potential source of leakage. 
Special measures should be taken to avoid all these defects. 
Traditionally mild steel reinforcement is used for water tanks 
because of the low allowable stresses permitted in the design 
to avoid cracks. But the superior bond properties of deformed 
bars, allowing for shorter laps and providing better crack 
distribution, will compensate for their sli^t extra cost even 
if their higher strength cannot be used to the fullest extent 
in the reservoir portion. However, the use of deformed bars 
will result in substantial saving in reinforcement required 
for the supporting shaft. Therefore, high yield strength 
deformed bars of 415 grade having characteristic strength of 
415 MPa have been adopted for the present work, 

3,4,2 Cracking phenomenon and its control 

In liquid retaining structures, serviceability limit 
state of cracking is the one which, in most cases, goveims 
the design of the reservoir portion and requires due consi- 
deration for crack ccaitrol in the design procedure. The 
complex and semi-randcxn phencxnenon of cracking has been 
tackled throu^ experimental work. Base et al, (1966), 

Beeby (1971, 197.9, 1983), Broms (1964, 1965), Broras and Lutz 
(1965), Desayl (1976a, 1976b, 1980), Nawy (1970, 1972) and 
Nawy and Blair (1971) are among the investigators who have 



contributed significantly to the present state- of -art. 

Although many expressions for estimating the flexiural crack 
width and crack spacing have been proposed , the significance 
of some of the basic variables is still being debated. There- 
fore, the procedure for crack control should be used only as 
a guide line rather than as a rigid rule. Formulae for flex- 
ural surface crack width are available in some of the latest 
codes. However, an equation which can be used with the confi- 
dence required for the design of water retaining structures, 
particularly when subjected to direct tension, is not yet 
available. Therefore, safety in respect of cracking due to 
direct tension, in water tanks, is ensutred only by designing 
the crack resistant structure by limiting the direct and 
bending tensile stresses both in concrete and steel to a 
certain level. To deal with such problems, BS 5337 specifies 
an alternative method of design in its clauses 6, 9 and 12, 

A similar design procedure, given in IS; 3370 (Part II) is 
adopted in the current study with a little modification in 
view of some latest provision in BS 5337. This procedure 
adopts the following measures to ensure waterti<^tness ; 

(i) Ihickness of the tank wall is chosen such that the 
direct tensile stress as well as flexural tensile 
stress in concrete, computed on the equivalent 
uncracked section basis, does not exceed the permi- 
ssible values 1.2 MPa and 1.7 MPa respectively for 
M20 concrete. 

( 11 ) The hoop reinforcement is designed to resist the 

complete hoop tension by restricting the pexrmissible 


stress in deformed bars to 100 MPa in accordance with 
•deemed to satisfy' provision of BS 5337 for almost 
no crack condition. 

(iii) Ihe reinforcement reqtiired to resist bending moment 
is based on the cracked section with the same value 
of permissible stress as given above in (ii). 

Although supporting shaft is subjected to combined 
axial load and bending moment and there is hardly any chance 
of crack width being critical/ it may sometimes require crack 
width calculations for dead weight + wind/seismic load xander 
empty tank condition. IS: 456 has given a guide line that the 
width of sxirface flexural crackS/ in general, should not 
exceed 0.3 rnm for mild exposure. CP 110 has given the foll- 
owing equation for estimating the width of surface flexviral 
cracks vhich is used, in the present work: 



(3.22) 


is the average strain in the member at the point of surface 
under consideration and is given by equation.: 


e 


m 



1.2 b^hy 



X 10 


(3.23) 


where a _ ** distance from the point xmder consideration 

to the surface of the nearest longitudinal 
bar, 

a= minimum cover to the tension reinforc^nenfc, 
h s depth of the monber , 



djj = depth of neutral axiS/ 

Cg = strain in steel at service loads,, 

y * distance from neutral axis to level of point 

considered, 

d = effective depth of the member, 

= width of the member at centroid of tension 
reinforcement, and 

fy = characteristic strength of reinforcement. 

The second term in Eq, (3,23) represents the stiff- 
ening effect of concrete in the tension zone. The allowable 
value of surface flexural crack width, for shaft, is taken as 
0,3 mm in the present study. 

3.4,3 Design of top dome, inner cylinder and ring beams 

It is found from actual analysis that the membrane 
and the bending forces in the top spherical dome and the inner 
cylinder are of such order that they do not demand thick 
secticxis. However, practical considerations preclude consi- 
deration of any thickness less than 75 mm for the top dome 
and 125 ram for the inner cylinder. The area of reinforcement 
is also governed by the minimum value of 0.25% except near 
the edges, in seme cases. The edge forces die out fast and 
extra reinforcement required at edges, if any, can be curt- 
ailed within a short distance. For economy, the rise of 
dome should be in the range of 0.25 to 0,4 times its radius 
and in order to avoid necessity of formwork for the top 
surface, the semi-central angle should not exceed 40® . Based 
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on these considerations the rise of the dome is kept as one- 
third of its radiiis. Cost of the dome and the inner cylinder 
in the objective fxjnction is computed, for a 75 ram thick dome 
with rise equal to one- third of its radius and 125 mm. thick 
cylinder and the reinforcement of 0.25% in both directions 
for these two elements except for some extra reinforcement in 
cylirxirical portion for a short length near junction 2, in 
some cases. 

Althou^fi; it is possible to provide the required hori- 
zontal reactions at junctions without providing ring beams, 
such an arrangemofit results in high tension in the edge regions 
and often requires heavy reinforcement and increased thick- 
nesses of the elements meeting at these junctions. Moreover, 
in liquid retaining stjructures, it is not desirable to have 
regions of high hoop tension, particularly at junctions, as 
theif may give rise to wide cracks. Therefore, ring beams 
provided at junctions 1 and 2 (see Figure 3,3) eliminate all 
such undesirable effects and can be easily reinforced. 
Dimensions of ring beams can be chosen as design variables, 
but cost of ring beams in relation to the cost of the system 
is very small. Besides efficiency of any optimization tech- 
nique will decrease with the increase in the number of design 
variables. As such, both the ring beams are chosen to be 
square in section with 1.5 t^ as the side, t^ denoting the 
thickness of the conical reservoir. Reinforcement required 
for ring beams is obtained from the hoop tension consideration. 



3.4.4 Distribution and arrangement of reinforcement in 
various elements 

Usually distribution and arrangement of reinforcement 
depend upon the forces in the structtire at service loads. As 
discussed in the preceding section the top dome and the inner 
cylindrical portion will have a xiniform distribution of 
reinforcement in both the directions except near the edges 
vjhere some extra reinforcement may be required, in a few cases. 

In order to control the cracks in the conical reservoir 
portion at service loads more hoop reinforcement should be 
placed at locations where hoop tensions are larger as obtained 
from the elastic analysis. The Figure 3.14 shows the actual 
general pattern of variations in the values of stress resul- 
tants- along the wall of the COTiical tank in the x-direction 
for reservoirs of different capacities. It can be seen that 
the hoop tension is larger in about two- thirds length of the 
lower portion except for a very short loigth near the narrow 
«id of the cone. Thus, the curtailment in hoop reinforcement 
for a short length, that too, at the narrow end of the cone 
will not be of much use. Whereas, in about one- third laigth 
of the upper portion of the cone, hoop tension decreases 
rapidly. Therefore, curtailment of reinforcement in this 
portion is meaningful and this results in considerable saving. 
The cottmon practice of curtailing the reinforcement where 
forces are smaller. Leads to the stepped distribution. 

Although the geometric parameters fixing the positions of 
points of curtailment can be taken as design variables, but 
this will increase !f.he size of optimization problem considerably 
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without much gain than those fixed by intxaitive trial and 
error method, it is important to note here that hoop tension 
variation for different tank capacities is found to be very 
much sJjnilar to that given in the Figure 3.14(a) and for this 
variation the distribution of hoop reinforcement given in the 
Figure 3.15/ which also indicates the points of curtailment/, 
is found to be the eccaiomical one for most of the cases and 
hence is adopted throughout the present study. The percentage 
of hoop reinforcement in the different regions and the thick- 
ness of the cone are taken as design variables in the optimi- 
zation process. The membrane and bending forces/ in the 
conical portion in x-direction do not demand the reinforcement 
more than the minimum of 0.25% except near the edges, in some 
cases. HerKie, in x— direction, the conical portion will have 
a uniformly distributed reinforcement throughout except at the 
edges where some extra reinforcement may have to be provided 
in scxne cases. 

For the supporting shaft, the area of reinforcement 
required in vertical direction will depend upon the holding 
moment and the axial force acting on any horizontal section 
of the shaft. These forces will be maximum at the section 
just above the foundation and minimum at the top of the 
shaft. A stepped distribution of reinforcement with points 
of eurtailraents at one-third and two-thirds heists of the 
ihaft as given in the Figure 3.16 is adopted throughout the 
current study* The percentage of reinforcement in vertical 
direction in three regi<xis and the thickness of the shaft 




Fig. 3*15 Distribution of hoop Fig. 3*16 Distribution of vertic 

reinforcement for the reinforcement for the 

conical tank * supporting shaft 
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are taken as design variables. A ininimum reinforconent of 
0.25% is fomd to be adequate in circumferential direction. 

In conical# inner cylindrical and supporting cylin- 
drical shells reinforcements in two directions are placed on 
both the sides# 50% on each face# in a staggered fashion. 

This is to take care of the sudden jumps in the bending moment 
(from positive to negative) near the junction 2, Any imper- 
fection due to construction by way of hump or necking will be 
taken care of by providing reinforcement on both the faces. 

The cover provided for conical and the inner cylindrical 
shells is 40 mm and for the support it is kept equal to 30 mm, 
the hoop reinforcement is placed on the inner side of the 
meridional reinforcement. 

3,4.5 Openings 

3. 4, 5,1 Stress concentration 

Small and big openings are required to be provided 
in water towers for different purposes. These openings res\ilt 
in weakcaning of the strxactural element locally at the points 
of their Icacation and in stress concentration around the 
edges of the opening. The latter is more serious foi^ small 
openings and the value of stress concentration factor (SCF) 
depmida upon the shape of the opwiing and the general stress 
pattern occurring in the structural elemait at that location. 
For uniaxial stress pattern of xaiiform intensity# in case of 
a small circular opening# the maximum stress is as high as 
three times the value of the average stress and it occurs at 



the ends of the diameter of the hole which is at ri^t angles 
to the uniaxial stress direction. Whereas, at the ends of 
the diameter parallel to the iiniaxial stress direction, a 
stress equal to the average stress but of opposite nature 
occurs. These stresses damp out rapidly within a distance 
equal to the radius of the opening and become equal to the 
average stress at a distance five times the radius of the 
opening as shown in the Figure 3.17. For biaxial stress 
pattern the stress concentratiai factor can be obtained by 
superposing the foregoing results of the uniaxial case. In 
large openings such as the one for the door provided at the 
bottom of the supporting idiaft for access to the reservoir 
portioi, the stress concaatration. factor at mid points of 
vertical edges of the opening can be obtained by approximating 
the opening as an elliptical one for which SCF = 1 + 2 ~ , 
where 'b' is the width of the opening and 'a* is its height. 
For an average size of door opaiing the value of SCF will be 
between 1.5 to 2 for uniaxial stress pattern of uniform inten- 
sity acting parallel to the axis of the shaft, whereas, the 
stress at the mid points of top and bottom edges of the 
opening will be equal to the intensity of uniaxial stress but 
of opposite nature. 

3, 4.5, 2 Reinforcing of openings 

A sufficient material should be added to compensate 
the weakming effect of the opening as well as to take care 
of stress concentration around the edges of the opening. The 
reinforcing material he placed immediately adjacent to the 




opening but suitably disposed in a profile and contour so as 
not to introduce an over-riding stress concentration itself. 
Reinforcement of an opening cannot be unquestionably obtained 
by adding huge amounts of material because this has the reverse 
effect and creates a ' hard spot' on the shell which does not 
allow natural growth under general strain pattern occurring 
througiKJut the shell elanent. This results in local overstre— 
ssing causing the failure of the structure which is similar 
to the pinching of a baloon. UsToally^ in reinforced concrete 
structures, best results are obtained by gradual increase in 
thickness of the concrete section and in reinforcing steel. 

In case of door opening in the vertical supporting shaft of 
a water tower the horizontal edges should be thickened with 
gradual increase in thickness so as to provide lintel action 
to take care of vertical forces at each of the edges. These 
edges are recjuired to take a small amount of vertical load as 
most of it will be transferred through the vertical reinforce- 
ment which will be bent to pass by the sides of the vertical 
edges of the opening. The amount of the reinforcement along 
the horizontal edges should be increased to take care of the 
tensile stress produced due to stress concentraticai effects. 

The veartical edges of the opening are stiffened by thickening 
the edges with ^adual increase in thickness so that the 
values of effective area of cross section and moment of 
inertia become at least equal to the respective values of a 
horizontal cross section of the jdiaft without opening. As 
far as possible, sharp corners should be avoided and the 



opening can be made approximately elliptical in shape so as to 
avoid stress concentration at comers as well as to facilitate 
the gradual bending of vertical bars to pass by the sides of 
the opening. 

3 . S Limit Analysis 

Although/ the limit state of collapse/ in this parti- 
cular case/ will not usually govern the design, but in order 
to find out the margin of safety against collapse, the limit 
analysis is carried out. 

3,5,1 Applicability of limit analysis 

The general theory of limit analysis of rigid-plastic 
structures has been rigorously formulated by Prager (1952). 

The theory of limit analysis can as well be appjlied to elastic- 
plastic or brittle-plastic materials. This woxild be possible, 
if it is first established experimentally that sucdh. structures 
deform considerably at practically constant load intensity 
before collapse. It shoxild also be possible to formulate the 
conditions of failure. Under-reinforced structures fall under 
this category. In the case of reinforced concrete shells, 
because the percentage of reinforcement is small, plastic 
b^avioir is fairly well established. Experimental investi- 
gations, by Baker (1953), Ernest and Marl ette (1954), 

Sawczuk (1961) and Das (1983), strongly indicated that 
collapse mechanisms and limit loads do exist for reinforced 
concrete shells. Adidam and Subramanyam> (1982) applied limit 
analysis to obtain collapse water pressure for reinforced 


concrete cylindrical tanks. They used lower boxind approach 
adopting the failure criterion based on the work of Sawczuk 
and Olszak (1961), whereas Adidam and Kalwar (1984) obtained 
the collapse water pressure through upper bound approach for 
a case of reinforced concrete conical water tank. This method 
is used in the current study also for finding the collapse 
water pressiure of the conical tank. 

Water tower considered in this chapter is mainly 
composed of four elements, viz., top spherical dome, conical 
tank, inner cylindrical shell and the supporting shaft. 
Although, the limit analysis for each element can be carried 
out, the possibility of collapse of the top dome or the 

inner cylindrical portion prior to the collapse of the conical 
tank or the supporting ^aft is negligible. This is because 
of the fact that both top dome and the inner cylindrical 
portion are primarily subjected to membrane compression and 
will collapse at a higher load factor than for the other two 
elements. Moreover, the possibility of combined mechaniam is 
also rare. Therefore, in the present study, limit analysis 
of ccatjical tank for axisymmetric case of loading that is dead 
load + water pressure under tank full condition is carried 
out. The supporting shaft can easily be analysed as a canti- 
lever subjected to combined axial and lateral forces. 

3.5.2 Limit analysis of conical tank 

As already mentioned, the limit analysis of the 
conical tank is carried out as per upper bound approach. The 
water pressure at collapse is obtained by considering all the 



possible collapse mechanisms of the tank. The design para- 
meters of the tank are fixed up by the serviceability require- 
ments. The Figure 3,18 shows the details of the conical tank 
with hydrostatic pressure loading at collapse and the ultimate 
hoop tensicxi capacity distribution based on serviceability 
limit state of cracking (see Section 3,4.4). The relative 
values of ultimate hoop tension capacities 
Ij^ arxi 1^, will depend on the ^ell parameters, 

3 . 5 . 2, 1 Collapse mechanism 

In this case, the conical reservoir is provided with 
a roof and a supporting structure, the ring beams and other 
el«nents meeting at junctions usiaally will not allow the 
conical portion to expand or contract considerably at the 
junctions and the only possible collapse mechanism will be of 
the form as shown in the Figure 3,19, The hinge circle 1 
will be always just near the bottom junction and obviously 
remains in the region ) . The hinge 2 will also be in 
regicm but hinge circle 3 may be formed in one of the 

three regions depending upon the shell parameters. Based 
upon the position of hinge circle 3, three different expre- 
ssions for internal work will be obtained. However, the 
expression for external work, due to hydrostatic pressure 
and self weight of the shell, will essentially remain same for 


all the cases 
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3 . 5 . 2 . 2 Internal work 

The total internal work is comprised of the energy 
absorbed at plastic hinge circles and in the meridional yield 
lines / energy dissipated by extension of the circumferential 
reinforcement. In a conical reservoir the value of meridional 
thrust is considerable and in some cases it may reduce 

the ultimate moment of resistance of the shell section. In the 
Figure 3,19, Mq 2 ^03 appropriate ultimate 

mcxnents of resistance at the hinge circles 1, 2 and 3 resped- 
tively. The expression for the internal energy absorbed at 
plastic hinge circles can be written as under 

^wl " "ox a ' i "02 (c *^03 j °°° “ 

The internal work by extensicwi of circumferential 
reinforcement will deporid upon the position of hinge circle 3. 
In the Figure 3.19, O is the origin, distances and X 2 
define the position of bottom and top of the conical shell, 

6^ is the displacement of any point X on the shell surface, R^, 
R 2 and R^ are the radii at top, bottom and point X of the 
conical shell. The point X will occupy the new position X' 
after virtual displacement and thus the increase in the 
length of the radius at X can be given for the two regions as 

6^ * 6^ sin a = — — ^ sin a f or x, < x < x- + a 

fOC. X ■ ■ ■ ^ <mm mm 2, ■ 


and 


(3.25) 





X. -f C - X 

* sia a ss — — sin a 

X c — a 


for x^+a£X£x^+c (3, 26) 


llius t±ie correspcwiding increase in circumferaitial length at 
this point can be given as 
2ii: (x - X.) 

6i » sin a for x, < x < x- + a (3,27) 

a 1 — — 1 

and 

2Tt (x. + c - x) 

61 «B — — sin a for x- + a<x<x- + c 

(3.28) 

Using the appropriate value for 6l given by Eqs. (3.27) 
and (3,28)# the expressions for inteimal work due to extension 
of cireuBmf€a:ential reinforcement for the three cases# depen- 
ding upon the x^ositicn of hinge circle 3 (Figure 3,20), can 
be written as 


for I2 5 c < X2 

^1"** (x - X ) 

Iw2 -/ sil. 0 . d^c 

(x, + c - a) 

*S ic - ai sinadx 

xj+a 

^l'*'^2 (x- + c - a) 

(x, + c - a) 

f 2 7iN^3 


sin a dx 


(3.29) 




for ^ c <_ 


Xj^+a 


^w2 “ / 


(x - X, ) 


tl a 


sin « dx 


x^+l^ 


- - (x. + c - a) 

* - ' ic - a) 


x^4a 


X^4C 


/. -ferry 


(x, + c - a) 


Xi+lj 


sin a dx 


(3.30) 


for a < c < 1. 


(x - X ) 

\,2 ~ f 2" J— i- sin a dx 


(x. + c - a) 

On simplification, the Eqs, (3,29) through (3.31) 
will reduce to the following three equations. 


(3.31) 


Oase I (1^ ^ c ^2^ 

K 

^w2 " ^ 2 “ + S-I^ - ®>“tl 

+ (ij - ii)C2c - I 2 - ijiNta + (d - 12>^ 

Case II ili < c < ^2^ 

^«2 “ *^2 * ”tl ■" ■ “’"tl 

+ (c - Ij)^ N^2] 


(3.33) 



Case III (a < c < 1^) 


^w2 * ^2 c (3.34) 

where K 2 * % sin a 

3 , 5 , 2 , 3 External work 


In the Figure 3.18# is the maximum hydrostatic 
pressure intensity at collapse. The pressure intensity p^ 
at any point X can be obtained as 


ox 


X2 - X 
^2 " 


where 


?o ® ” ^ 1 ^ ' 


P »= unit weight of water# and 

= partial safety factor for loads. 

The expressicxi for the work done by hydrostatic pres- 
sure through virtual displacement (Figures 3,18 and 3,19) can 
be obtained as follows 


2 itp^ cos a .^l'** 


VI 


y C / x(x 2 - x)(x - Xj^) dxj 

2 ' X ' iX*! 


2 % cos a 


X^4C 


(S' - - X-'V [ •/' - ^>‘*1 + c - X) dx] 

^ i X« *4^ - 

(3.35) 

On simplification of Eq, (3,35) the follcwing expre- 
ssion is obtained 


\l “ '^3 ‘ 


2 2 

(X 2 + Xj^Xa + 3Xj^ + 3Xj a) 
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3 2 2 3 

a + 4a + 6ax^ + 4x^ x^x^Ca + 2Xj^) 


^ ^2^1 ^^2 ^ ^ ^ + a + ) ^2 ^ 


tc + a^ + ca + 3Xj^(c + a + x^)] + -j [(c + a)(c^ + a^) 

(3, 36 ) 


+ 4Xj^(c^ + a^ + ca) + 6xJ(c + a) + 4xJ] } 


where K, 


2iip^ cos a 
^2 ~ ^1 


If is the \in±t weight of concrete/ then the 
dead load at collapse w^ per \init surface area of the conical 
shell can be written as under 


w. 


Yr P t 
' L C C 


where t =» thickness of the conical shell, 

%# 

Proceeding in a similar manner/ external work due to 
self weight of the shell can be obtained. The final expression 
comes out as follows: 


®w2 • *^4 ^ 


2 2 

a + 3Xj^ + 3x^a Xj(a + 2Xj^) 


^"^1 ■*’ 4- a + 2x^) -K a^ -f ca 


- Xj^(c a + ] 


(3.37) 


where « 2 w^ cos a . 

ThuS/ the principle of virtual work gives 


^wl ^w2 


®wl ®w2 


(3.38) 



The Eq. (3.38) can be solved, for all the possible 
combinations of a and c, to obtain the value of p for each 
such combination. The lowest value of p^ so obtained will 
correspond to the hydrostatic pressure at collapse. In the 
present study one hundred and twenty one such mechanisms are 
tried to find out the collapse pressure p^. In Eq. (3.38), 
the value of (the external work due to the self weight of 
the shell) is £o\ind to be negligible as compared to that of 
®wl* contribution of to total external work is hardly 

5^. Similarly the contribution of I^^ (the internal energy 
absorbed at plastic hinge circles) found to be about 10% of 
the total internal energy. Hence, for simplificatican , values 
of E ^2 neglected in finding the collapse 

pressure p^* 

3,5.3 Limit analysis of cylindrical supporting shaft 

In case of small capacity tanks, the design of the 
shaft is foxind to be governed by serviceability limit state of 
deflection* Whereas in some cases of tanks having capacities 
500 kl to 1000 kl^ ultimate limit state of strength of a local 
section of the cylindrical Shaft when considering an element 
of unit width, may be critical. However, the limit state of 
collapse of the shaft is never found to govern the design. 

In order to estimate the margin of safety against collapse, 
the limit analysis is carried out by idealising the shaft as 
a cantilever of hollow circular section, it is of interest 
to note here that under combined axial and lateral forces, 
a local section of shell element of unit width will have 



almost uniform axial strain across the thickness of the shell. 
Therefore, to safeguard against failure of such local sections, 
the maximum strain in the extreme fibre of the hollow circular 
section should not exceed 0,002 as per IS: 456. 

3,5.4 Material properties and assumptions 

In order to find out the ultimate strength of a 
section xinder axial, bending or combined axial and bending 
forces, the stress-strain relations used, for both concrete 
and steel, are as given in IS; 456, For obtaining the design 
strengths o'f these materials the partial safety factors speci- 
fied by the Code are 1,5 and 1,15 for concrete and steel 
respectively , 

The following assumptions, as per IS: 456, have been 
emoloyed in the limit analysis: 

(a) Plane sections normal to the plane of bending remain 
plane after bending, 

f 

(b) The strain in concrete at the outermost compression 
fibre in flexure reaches a value of 0,0035 at failure, 

(c) The distributich of compressive stress in concrete is 
defined by an idealised stress-strain curve as shown 
in the Figure 3, 21(b), 

(d) The tensile straigth of concrete is ignored, 

(e) The stress in the reinforcement is derived from the 
represaatative stress-strain diagram of the steel, which 
is shown in the Figure 3, 21(a). 

(f) The roaxiJWM compressive strain in ccwicrete in axial 

compressicai is 0,002 at failure. 
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(g) The maximtjin strain at the most highly compressed 

extreme fibre under axial load and bending shall be 
taken as 0,0035 minus 0,75 times the stiain at the least 
coftitMTessed fibre when there is no tension on the section. 

3 . 6 Ootimization 

Most of the engineering designs turn out to be non- 
linear constrained optimization oroblems which can be stated 
in the standard mathematical form as: 

Find vector X of the design variables which minimises 
the objective function F(X) subject to m constraints 
gj(x) <0, j * 1, 2 , ... m. 

The objective function considered in the present study 
is to minimize the cost of the superstructure of the water 
tower. Ihe choice of the design variables arid the constraints 
on the design are discussed in the subsequ«it sections with 
special reference to certain important asnects of the problem 
at harxi. The -orecise mathematical programming formulation 
subsequently follows, 

3,6,1 Organisation of optimization process - some practical 
conslderat icai s 

Optimization process is a cyclic design-analysis 
procedure. The computational time to obtain the optimum 
design increases with the increase in the nrjmb^ of design 
variables and the nwber of constraints in the pxX^Irem. 
Therefore, a serious thought is required to be given to the 
choice of desi^ variables and the ccaastraints to be included 



in the optimization framework, with the variation in the 
design variables, the design forces also vary. However, if 
the change in the design forces is very small with regard to 
the variations in the design variables, the design forces can 
be considered to be invariants in the optitnizaticxi process. 

In that situation the design forces can be given as a slightly 
conservative inout and this avoids the computation of the 
design forces in each cycle of optimization. This leads to a 
substantial saving of computational time in the overall optimi- 
zation effort. Although such a situation is a rarity, but in 
the tsresent study it is so and this has been established while 
observing critically the results of the elastic analysis for 
finding the design forces in the reservoir iDortion, In view 
of the foregoing considerations, certain important aspects to 
be incorporated in the optimal design process are discussed 
herein . 

Ihe governing design forces in the reservoir portion 
are obtained from the elastic analysis for axisymmetric loading 
combination of dead load and water pressure loading xinder 
tank full condition. These forces are the functions of 
material properties and design parameters like , radius of 
the supporting shaft '^2 * angle 'a' and the thicknesses 
of the various elements of the tower. The influence of 
these variables cffi the design forces of the reservoir portion 
is found to be as under: 

(a) Figures 3.22 and 3.23 show the variation of the stress 
resultants in the reservoir portion with a and R 2 res- 
pectively, for a typical bank of 300 kl capacity, it 
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.3-23 Variation of stress resultants along the wall of the coni 
tank with R? 





can be seen that the forces change appreciably with 
change in both the parameters . 

(b) For a particular value of R 2 and a there is hardly any 
change in the governing design forces, particxalarly 
that of the hoop tensions, with the change in thicknesses 
of various elements. Although the bending moment 'M ' 

Jv. 

changes slightly which usually does not govern the 
design, it is limited to edge region only. The Table 
3 , 17 indicates the variation in the governing design 
forces with thickness of conical tank and the shaft, 
for a reservoir of capacity 300 kl . 

Unlike the reservoir portion, the governing design 
forces of the supporting shaft vary appreciably not only with 
Rj and a but also with the change in the thickness of the 
shaft. 

In the light of the above observations, there can be 
two possible approaches to achieve the optimal design. One is 
the direct approach and the other is an indirect one. 

In the direct approach R 2 and a could be considered 
as design variables with other possible variables like thick- 
nesses and reinforcements for different elements of the tower. 
The optimal design so obtained will directly give the optimal 
shape parameters , R 2 and a, and corresponding values of the 
other design variables. 

In the indirect approach R 2 and a can be made as 
jareassigned parameters and the optimal design so obtained 
will give the values of corresponding design variables for a 
particular combination of and «, but with a fixed input of 



'I^ble 3,17 Variation in governing design forces of the reservoir portion with thickness of 
conical tank and the aupptJrting shaft 
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design forces for the reservoir portion which are slightly 
on the conservative side. The design forces corresponding to 
the optimal design so obtained can finally be evaluated for 
the reservoir portion and the optimal design rechecked for 
the actual design forces. The difference between the actual 
design forces and the input design forces is observed to be 
negligible basically because the membrane forces which govern 
the design are insensitive to the thickness variations. There- 
fore# the revision of the optimal design so obtained is never 
necessitated. In order to obtain the optimal configuration 
of the tower through indirect approach# a parametric 
study for different combinations of R 2 and a is made in the 
present work. In order to arrive at the optimal configuration 
atleast nine different sets of R 2 and combination have been 
xised in the jaresent work in each case. 

The Table 3.18 shows the results of optimal designs 
obtained through the two approaches# described in the fore- 
going discussion# for the following design data? 

3 

Capacity of the tank - 100 m 

Height of the staging *= 20 m 

Tree board ** 150 mm 

2 

Wind 7x>ne » II (Basic wind pressure = 1.5 kN/m ) 

It can be seai from the results given in Table 3.18 
that the optimal designs obtained by the two approaches are 
almost the same but the total computational time required 
throu^ indirect approa<^ is very much less than that required 
by direct approach. Moreover# the indirect approach gives 



Table 3,18 Results of optimal design through two approaches 
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optimal designs for different ccxnbinations of R 2 and a in 
addition to the one which refers to the optimal configuration. 
As mentioned earlier , in the present study , the objective 
function does not include the cost of the foundation, there- 
fore, the optimal diameter so obtained by direct approach may 
not sometimes suit the foundation conditions for overall 
economy. In such cases one has the choice to go for an 
optimal design with a suitable diameter of the supporting 
shaft, out of those which are obtained in the process of 
indirect approach. Thus, indirect approach has many advantages 
and hence is adopted herein. To arrive at the optimal confi- 
guration through indirect approach values of and R^/R^ are 
increased or decreased by 0.25 m and 0.5 respectively. These 
intervals are found to be satisfactory both for arriving at 
near exact optimal design and for standardisation of such 
designs and their formwork to be used. 

It is not always necessary to include all the design 
requirements as constraints in the optimization process. In 
order to save conputational time, only those constraints 
should be included which are likely to govern the design. The 
optimal design so obtained can finally be checked for other 
design requirements and may be modified, if needed. In view 
of this consideration, the constraints not included in the 
optimization process, but almost invariably found to satisfy 
the final optimal design, are as follows: 

(i) The constraints related to safety against collapse 

for reservoir portion as well as the supporting shaft. 
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(ii) The constraints dealing with limiting surface flexural 
crack widths in the shaft. 

(iii) Some of the constraints which deal with minimum require- 
ment of reinforcement as eiqjressed by the inequalities 
(3.40)# (3.41) and (3.47)# which are given in the 
, Section 3, 6 , 4 . 

It does not seem feasible to carry out optimal design 

study# in the present case# for all possible combinations of a 

wind zones and seismic zones defined in the country without 

prohibitive expenditure of computer time. Moreover# it has 

been recognised that designs for all such combinations are not 

indeed required to be carried out because a water tower 

designed for a particxAlar wind zone may be safe in many 

seismic zones. For example# a tower of 100 kl capacity and 

20 m staging height# designed for wind zone I (basic wind 

2 

pressure 1 kN/m ) is found to be safe for all the five 
seismic zones of India. Therefore, in the present work# 
optimal designs, are carried out for the various wind zones 
of India and their sxaitability is established for seismic 
zones. Thus# the analysis for finding the seismic forces is 
carried out only at the time of checking the suitability of 
the final design for a particular seismic zone. This results 
in considerable saving of .ccxnputaticaial time. 

To stmmarise the above discussion# the optimal design 
of the siiperstructure of the water tower# in the present work# 
is carried out for different wind zones by performing a 
parametric study for a set of values of R 2 and oc. The optimal 
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design so obtained is finally checked for its safety and suit- 
ability in various seismic zones. 


3.6,2 Design variables 


The vector X gives the various design variables consi- 
dered in the present design problem, which is expressed as 
under : 


^rcl 

Prc2 

Prc3 

Prc4 

^s 

^rsl 

^rs2 

^rs3 

Prey 


where t ^ thickness of the conical tank wall, 

c 

Pj.^ 2 Prc3 percentages of hoop 

reinforcemaat in the three regions of the conical 
tank (Figure 3.15)/ 

p , s= percaitage of reinforcement in meridional 

direction for conical tank, for a short length 
at junction 2, 
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tg = thickness of the supporting shafts 

^rsl' ^rs2 ^rs3 percentages of reinforce- 

ment in vertical direction in the three regions of 
the supporting shaft (Figure 3.16)^ and 
Pj^cy “ percentage of reinforcement in x-direction for 

inner cylinder, for a short length at jxmction 2, 

The values of maximum bending moments, for ccxiical 
and inner cylindrical portion of the reservoir, occur in the 
edge region of jianction 2, However, these moments reduce at 
an expcxxential rate and hence reinforcement p . and p need 

to be provided only for a short length near the junction. For 
the remaining length of both the elements, a minimum reinfor- 
cement of 0.25% has been found to be more than adequate. Any 
of these lengths is not forond to be more than 1 m, hence in 
computing the objective functltxi the reinforcanents Pj-j,^ and 
p are ccxisidered for 1 m length only in the two elements. 

3,6.3 Objective fmction 

The cost of concrete, reinforcanent, formwork required 
for the water tower including labour but excluding the found- 
ation, is dhosan as the objective function, in the present work. 
Adopting the design for top dome and inner cylindrical 
portion as givai in Section 3.4,3 and the reinforcement distri- 
bution for conical tank and support as mentioned in Section 
3.4.4, the expression for objective function can be written as: 
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F 


> 4 iere 


2 7 rR^<cog^ « )(t,c + ^ ^ 0 ^25 

X a c in>r“ 

1 *'2 


d'-c 

cosa~^ (tc + 2^ t 


^ (R? - R^) 


100 ^s ^s ■*■ 


■c c i6o~ c ) 


" (R? - R?) 


2' Pr 


,-,2 , 


* — ^ C^t C n , , " <R 4 - R3) ,R 


COS a 


100^ ^c ^S + (^^£2 . ^ ^ 

^ S cos a ^ 100 ^c ^' 


7E (r2 « 

+ - 1 ,Pro3 , 

cos a 'loo P c ) + 2 ir L (t c 

cy cy cy c 


s s 


+ ^ X 0. 2!^ 


100 ‘cy 's =3 2C ) + 27. R L [t c + 2i25 t „ . 

2 *-3 0 100 s ^ 


s s 


^<s C 

+ 2C 4- - ” 3 S / 

f 300 <Prsi Prs2 ■" P,-»3) ] 


-rs2 - Prs3 

^ C(R/-i 4* l^n r»oQ iv \ 2 v^2 


2J^1.0 003 0 .)"- R^] . 0.25 

^ ioo Po c^) 


cos a 


G s S' 


+ 2TtR r ^cy " 

cy^ 100 ^cy Ps 


<3,39) 


= radius of the top dome, 

R^y = radius Of Inner cylinder (which Is equal to R 
in this case), 

Cj, » cost of finished concrete In rupees per m^, 

C3 = cost Of steel In rupees per )<n, 

Cf = cost Of fomwork In rupees per square meter, 

‘cy = thickness of the top dome and the Inner 

cylinder respectively 

% » weight of steel in kN/m^, 

L^y and L . length of Inner cylinder and supporUng 

*aft respectively, and all other terms are as 
defined earlier. 
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3.6.4 The constraint s 


To take care of the temperature variations and shrin- 
kage effects^ the percentage of reinforcement for deformed 
bar*/ in each of the two directions should not be less than 
^r min' ^3 equal to 0.25% as per IS: 3370. ThuS/ the 

eight constraints can be expressed in the form 


and 


min 

^rcl 


- 1 < 0 / 


^r min 
Prc2 


- 1 < 0 


^r min 


■rc3 


- 1 < 0 / 


r min 
^rc4 


- 1 < 0 / 


^r min 

p 

'^rcy 


- 1 < O / 


■r min 


P. 


- 1 < 0 / 


rsl 


min 

Prs2 


- 1 < 0 / 


r min 
^rs3 


- 1 < 0 . 


( 3 . 40 ) 

( 3 . 41 ) 

( 3 . 42 ) 

( 3 . 43 ) 

( 3 . 44 ) 

( 3 . 45 ) 

( 3 . 46 ) 

( 3 . 47 ) 


Practical ccxisiderations dictate a minimum thickness , 
for conical tank wall and the cylindrical shaft, 'mere 
is generally no advantage in adopting a thickness less than 
1SK3 nTO-125 ram as the precision required in fixing the reinfor- 
cement/ placing and compacting coiwrete would involve extra 
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cost. An error in the displacement of the reinf orc^noat/ 
provided to resist bending moment/ will be more serious in 
thinner walls, A minimxjm thickness of 125 mm is adopted for 
both the elements in the design of Type I towers considered in 
this chapter. Therefore/ the ninth and tenth constraints are 
written in the form 


min 


- 1 < 0 / 


(3.48) 


and 


min 


- 1 < 0 . 


(3.49) 


At service loads cracks due to hoop tension in the 
conical reservoir portion are controlled by limiting the 
direct tensile stresses in steel as well as concrete to the 
maximum permissible values/ f and f '^l^ch are 

100 MPa and 1,2 MPa for steel and M20 concrete respectively. 
These Checks are to be carried out for each of the three 
regions of hoop reinforc«nent. ThuS/ the constraints eleven 
through sixteen cab be written as 


(Nqi) 


max 




St max 


] - 1 


(3.50) 




'#2' ' '^c^ lOO St max 

max 


(3.51) 


^^03 V /^^c^lOO ^ ^st max^ “ ^ 1 ° ' 

max 


(3. 52) 


^^01^ 


max 


__ 
eul ct max 


- 1 < 0 


(3,53) 
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A _ 
eu2 


I 


) 

max 
ct max 


1 < 0 , 


(3.54) 


and ^ 1 < 0 , (3.55) 

eu3 ct max 


^Qn± equivalent uncracked sectional area of concrete 

in ith region of hoop reinforcement which is given by 


eui 


s= t. 


^ (m ~ 1)P--. 

fl + —.■■I . I ■— . i £S:4: , I 

100 -J 


Where m « modular ratio. 


Serviceability limit state of cracking also demands that 
flexural tensile stress in concrete^ computed on equivalent 
uncracked section basis, and in steel based on cracked 
section ar# less than 1.70 MPa and 100.0 MPa respectively. 
Maximvim bending moments in x-direction for both conical and 
inner cylindrical shells occur near the edge region at juncticxi 
2. However, in ccxiical shell the value of stress resultant 
N will also be maximum at juncticxa 2 which will nullify the 

.wV 

bending tensile stress and there will be hardly any need for 

carrying out the check for flexural cracks. Whereas in case 

of inner cylinder N is caused only due to self weight of the 

cylinder and for all practical purposes it is neglected. If 

M and M are the maximum bending moment carrying capaci- 
cyc cys 

ties of the section of unit width for the inner cylinder, as 
per the limiting stresses defined in the foregoing for 
concrete and steel respectively. Then the constraints seven- 
teen and eighteen can be expressed in the form 
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M 

- 1 

M ^ 

eye 


(3.56) 

and 

M 

.91^ . 1 

M 

cys 

• 

o 

V| 

(3.57) 

where 

^cya ~ actual maximum bending moment in the edge 



region at junction 2 for ijuiit width along 



circtomference 

of the inner cylinder, at 



service loads 

• 



nie limit analysis, carried out to assess the strength 
of the tank to resist the water pressure, furnishes the value 
of collapse pressure p^. Using a partial safety factor 1,5 
against collapse, the nineteenth constraint can be written as 


1.5 P {x .2 - sina 


- 1 


(3.58) 


If M. is the ultimate moment of resistance of the 

urs 

circular shaft at its critical section under combined axial 
and lateral loads at collapse and M is the actual value of 
the balding moment at that section, then the safety against 


collapse demands that 

^ urs — uas 

constraint may be stated as 


Thus, the twentieth 


M. 


vias 


M 


- 1 


(3.59) 


urs 


Although the collapse in real saise will occur when 
maximum water pressxire in the conical tank reaches the value 

p or the maxiimm bending moment at the critical section of 

O ' 
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the shaft becomes equal to M > but the structure will become 

vlXTS 

^serviceable as soon as the first hinge circle foms in the 
conical portion or considering an element of unit width# the 
xiltimate limit state of strength of a local section of the 
shaft is reached. Usually the limit state of collapse will 
not govern the design of the tower. Whereas the other two 
cases of failure# due to ultimate limit state of strength being 
reached at local sections# may be critical in some cases. 

Hence# safety against such failures is more important. The 
first hinge circle in the conical tank will be formed in the 
edge region of junction 2 because both M and N are higher in 
this region. Whereas the critical sections for the shaft will 
be one just above the foundation and other two at points of 
curtailment of vertical reinforcement. If is the ultimate 

moment of resistance per xinit width of the conical shell at 
the critical section in presance of axial thrust and is 

the actual value of maximum bending raomont at that section# 

thon for safety M > M, Thus the twentyfirst constraint 

can be written as 


M, 

M 


uac 

■\irc 


- 1 < 0 


(3.60) 


Under dead load + water load combination# when tank 
is full# the shaft is subjected to axial compressive stresses. 
For other loading combinations also# due to combined action 
of axial and lateral forces# a local section of supporting 
cylindrical shaft of unit width will have almost uniform 
axial strain across the thickness of the shell. Thus# the 
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ultimate strength in axial compression^ P ^ /> and 

ULITSl iUTSZ 

^urs3 such local sections is required to be greater than 

the actual maximum ultimate forces/ P , « P « and P 

uasl' \aas2 uas3' 

coming on to these sections,. Hence the constraints twenty-two 
through tw«ity-four can be written as 


and 


P 

uasl 

p 

ursl 


< 0 / 


P 

ua3 2 

P 

urs2 



p 

uas 3 

P 

urs3 



(3.61) 


(3.62) 


(3,63) 


The constraints twenty-five through twenty-seven 
dealing with surface flexural crack widths/ 

w^r 3 ^ due to maxinwm bending moments in the three regions of 
vertical reinforcement of the supporting shaft/ can be expre- 
ssed in the form 


w 


crl 


w 


- 1 < 0 / 


era 


(3.64) 


w 


cr2 


w 


- 1 < 0 / 


era 


(3.65) 


w - 

and rr^ - 1 < 0 . (3.66) 

era 

where w » the maximum allowable width of crack, which 
era 

is taken as 0.3 mm. 

As suggested in Sectiai 3.4.2/ the procedure for 
cracks control should be used as a guide line because its 
reliability is still being d^ ted. Therefore an indirect 
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way of crack control by lAmiting the maximum deflection still 
has its place in the design procedures. Moreover^ deflection 
of the supporting shaft# due to lateral loads# results in 
additional moments induced by vertical loads. If the deflec- 
tion is not restricted# these moments may considerably affect 
the design of the i^aft and the foundation and thereby the 
overall economy. Usually for tower- type engineering installa- 
tions the maximum allowable deflection is taken between "^q 
1 

to ISo height. However in the present study the 

maximum allowable deflection at the top of the shaft is taken 
as hei^t of the ^aft, Thus# the twenty -eighth 

constraint can be written as 


I. 


act „ 1 ^ 


(3.67) 


‘all 


where ^^ct ~ actual deflection at the top of the shaft# and 
^all ~ n^ximum allowable deflection. 

To safeguard against elastic instability of the supp- 
orting shaft, the procedure given in IS: 2210 is adopted in 
the present study. For safety the actual compressive stress, 
^act' service loads should not exceed the permissible 
critical stress, Thus# the twenty-ninth constraint 

is given as 


^22^ - 1 < 0 (3.68) 

*^crit “■ 

Although the amount of wind induced ring moment in 
the shaft is small and does not demand any constraint for 
safety and serviceability# but for completeness the thirtieth 
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constraint dealing with the surface flexural crack widths 
can be written as 

- 1 < 0 (3,69) 

'^cra 

3,6.5 Optimal design formulation and solution procedure 

From the foregoing discussion it can be concluded that 
the formulated optimal design problem to obtain the optimal 
configuration and the design of the water tower turns out to 
be a problem of ten design variables, as defined in Section 
3.6,2, and having thirty constraints as described in Section 
3,6.4. The objective function defined in the Section 3.6,3 
is a ncai-linear f\inction of the design variables and so also 
most of the constraint functions. Therefore, the formulated 
optimal design probl«n is a non-linear constrained programming 
problem v>4rLcai has been converted into unconstrained non- 
linear programming problem through the use of interior penalty 
function, Ihe unconstrained non-linear programming problem 
so obtained is finally solved as a sequence of unconstrained 
mljiimization problems by Davidon-Fletcher-Powell method as 
outlined in Section 1.4.5. 

3,7 Parametric Study 

Optimal limit state design of Type I water towers 
has been carried out for the following cases; 

(1) Capacity in m^s 100, 200 and 300. 

(2) Staging height in ms 15, 20 and 25. 

2 

(3) Basic wind pressure in MJ/m : 1.0, 1.5 and 2.0. 
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The suitability of these designs# obtained for various 
wind zones#^ is established for the different seismic zones, 

A free board of approximately 150 mm has been provided in all 
cases. 

3.8 Rates of Construction Materials 

Current rates of the different construction materials , 
which are used for the evaluation of the cost of the super- 
structure of the tower are given as under: 

3 

Cost of finished cement concrete of grade M20 = Rs, 750 per m 

Cost of high yield deformed bars including cost 

of cutting# bending and laying = Rs. 600 per kN 

2 

Cost of formwork = Rs, 600 per m 

In order to investigate the effect of escalation in 

the material cost# the practice of using cost ratio of roater- 

3 2 

ials is followed. Cost of one m finished concrete and one m 
of formwork have been normalized with cost of one of rein- 
forcemeant and the resulting cost ratios are termed CRl and CR2 
respectively, Thus# the values for CRl and CR2 as per current 
rates are 1,25 and 0,10 respectively. 

3.9 Results and Discussions 

The optimal design results have been presaited essen- 
tially in the tatJUlar form. For a given capacity of tank# 
staging height # seismic and wind zones an optimal design can 
be picked frcxn the design Tables 3,19 through 3.21 which 
gives tooth optimal configuration and design variables for the 


tower 
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Table 3.21 Optiiml designs of Type I towers with 300 m capacity tank 
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The Figure 3,22 shows the variation of stress resul- 
tants along the wall of the conical tank with a for the tank 
of 300 kl capacity which is supported on the shaft of 1.125 ra 
radius. It can be seen that for a = 0,957 radian the maximum 
hoop tension and bending moment are minimum and the shearing 
force also reduces considerably. The Table 3.21 gives 
optimal designs for 300 kl capacity tank. It is interesting 
to note that for some cases of wind zones and heights of 
shaft the optimal configuration also turns out to be a = 0,957 
radian and Rj = 1,125 m. Thus# the results of the Figure 3.22 
and Table 3,21 together indicate that the optimal configuration 
of the reservoir adjvists itself in a manner that the bending 
moment and shearing forces at junctions reduce considerably 
bringing the state of stress in the reservoir portiop in the 
close proximity to that of the membrane state. It is also 
observed that the value of maximum hoop tension in the reser- 
voir becomes minimxim for the optimal configuration# in most 
of the cases. 

The Table 3,25 indicates the effect of cost ratios on 
optimal configuration. It can be seen from the results that 
for the various corobinaticais of cost ratios CRl and CR2/ the 
optimal configuration comes out to be the same. It is also 
observed that for a particular shaft radius the optimal 
value of ^ 1/^2 same for different cost ratio 

coirfsinations . The Table 3.26 Shows the effect of cost ratio 
on optimal design variables. There is hardly any change in 
the design variables of the reservoir portion but the design 
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Table 3.23 OpUrral designs of Type I towers with 200 capacity tank for non-optimal shaft 
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Table 3.25 


Effect of cost ratios on optimal configuration 




Capacity of 

tank = 

300 m^, L = 

25 m, Pj^ = 1.5 

s. 

CRl 

CR2 

"^2 

Rt/R^ 

Optimal 

Optimal 

No, 





cost 

conf iguration 




(m) 


(Rs) 


1 

1.25 

0.10 

1.375 

4.5 

194073 






5.0 

190509 






5.5 

196318 

R^ = 1,625 m. 




1.625 

3.5 

192113 






4.0 

186975 

R-/R^ = 4.0 or 





4.5 

190509 





1.875 

3.0 

208645' 

a = 0.924 rad 





3.5 

194895 






4,0 

198848 


2 

1.T5 

0.10 

1.375 

4.5 

223060 






5.0 

218856 






5.5 

225521 

R« = 1.625 m. 




l4J^5 

3*5 

222096 






4.0 

211467 

R,/R„ = 4.0 or 





4.5 

218041 

X z 




1.875 

3.0 

237973 

a = 0.924 rad 





3.5 

222587 






4.0 

229977 


3 

1.25 

0.15 

1.375 

4.5 

227892 






5.0 

225197 






5.5 

233870 

^ l*62!iS 




1.625 

3.5 

228795 






4.0 

224301 

R-/R^ 4.0 or 





4.5 

230481 

X z 




1.875 

3.0 

248978 

a = 0*924 rad 





3.5 

235343 






4.0 

240764 
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Table 3.26 Effect of cost ratios on optimal design variables 

• 3 2 

Capacity of tank = 300 m ^ L = 25 Pj^ = 1.5 kN/m 


— 

r."'""’"" 

r- 








— n 


s. 

No. 


V«2 

CRl 

CR2 

Optimal design variables 

Optimal 

cost 






^c 

*^rcl 

^rc2 

r 

Prc3 

Prc4 








^rsl 

^rs2 

Prs3 

P 

rcy 



(m) 



^ - 


1- - - .. . . 

(mm) 





.. 

(Rs) 



1 

1.625 

4.0 

1.25 

0.10 

132 

1.435 

1.090 

0.676 

0.25 

186975 






143 

0.445 

1.233 

1.888 

0.25 


2 

1.625 

4.0 

1.75 

0.10 

131 

1.444 

1.080 

0.690 

0.25 

211467 






135 

0.610 

1.309 

1.865 

0.25 


3 

1.625 

4.0 

1.25 

0.15 

133 

1.390 

1.060 

0.695 

0.25 

224301 






143 

0.445 

1.235 

1.886 

0.25 


4 

1,625 

4.0 

0.75 

0.10 

133 

1,434 

1.047 

0.700 

0.25 

160955 






145 

0,429 

1.096 

1.858 

0,25 
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variables of the supporting shaft change slightly with the 
change in the cost of concrete. It is observed that with the 
increase in cost of concrete , for cost of formwork remaining 
the same, the thickness of the supporting shaft decreases aix5 
amount of the reinforcement increases slightly. Whereas for 
a particular cost ratio of concrete to steel the design vari- 
ables are not found to be sensitive to any change in the cost 
of formwork. Thus, for a specified cost ratio of concrete to 
steel, it is observed that both the optimal configuration and 
design variables are not sensitive to changes in the cost of 
formwork. On the other hand for various cost ratios of 
concrete to steel, only the optimal design variables of the 
supporting shaft do vary but only marginally, while the optimal 
configuration remains essaitially unaltered. 

It is clear from the foregoing discussion that changes 
in costs of different materials will not alter the optimal 
configxiration of the tower and the design variables of the 
reservoir portion and only the design variables of the supp- 
orting shaft change very slightly with the change in the value 
of CRl, Thus for all practical purposes the optimal design 
obtained using absolute costs of materials prevailing at the 
time of investigation will also be optimal or very near to 
the optimal design vdien the costs of materials are different. 

Ihe design Tables 3.19 through 3.21 show that in 
case of a particular capacity of tank, for more than caie 
wind zone and heights of staging, the optimal configuration 
comes out to be same. In such cases it is expected that the 
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design variables of the reservoir portion also remain the 
same as the governing design forces of the reservoir do not 
change with changes in wind forces and the staging heights. 
However, a little variation found in these variables is unavoi- 
dable in solutions based upon any nxmerical technique. Ihus, 
while standardising these designs, this particular fact is to 
be kept in view. 

In some of the cases, thicknesses of the conical tank 
and the supporting shaft have assumed the minimum value. In 
optimization studies, when such a situation is encountered, 
it is normal practice to explore vhether the specified minimum 
can be further reduced, in the present case, since a minimum 
nominal cover has to be provided and the value of 125 mm has 
been arrived after ccxisidering the limitations imposed by the 
construction aspects, it is not desirable to reduce the 
specified minimum thickness any further. 

irom design Tables 3.19 through 3,24 it can be 
cd^aerved that for small capacity tank the seismic forces are 
very small as compared to the forces due to wind and the value 
of seismic forces increases considerably with the increase in 
capacity. For 100 kl capacity tank the design obtained for 
wind zone I itself is very safe even for seismic zone v, for 
all staging heights. Whereas, in case of 300 kl capacity 
tank, seianic forces in zone iv are comparable to the wind 
forces of zone I for 15 m staging height. It is also obser- 
ved that for a given capacity of tank, wind and seismic zones, 
the ratio of the seismic forces to wind forces increases with 
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decrease in staging height, ThuS/ in some cases, for small 
staging height, the seismic forces may govern the design of 
the supporting structure whai the capacity is more. 

Ihe design Tables 3,22 through 3>24 give optimal 
designs for a radius of supporting shaft whi<dv is in the 
neighbourhood of the optimal radius. These designs are in 
addition to those corresponding to the optimal radius and may 
come out to be economical, for a particular foundation condi- 
tion, when the cost of foundation is also taken into account. 
It is interesting to note that the optimal value 
of oc lies between 0,779 to 0,977 radians for all the cases 
and between 0.85 to 0,96 radians for most of these cases. 

The forces in the ring beams at junctions 1 and 2 
are of such order that 200 x 200 mm ring beams with 0.50 per- 
cent hoop reinforcemait is found to be more than adequate for 
all the cases. 

Constraints which are critically satisfied are those 
regarding the following: 

(1) minimum thickness of tank wall, 

(2) minimum thickness of supporting shaft for some of the 
cases, 

(3) minimum valxie of reinforc«nent in x-direction of 
conical tank, 

(4) maxinwra values of stress in the hoop reinforcaments 
in all the three regions of the conical tank, 
maximum value of direct taasile stress in concrete in 
the region of roaxiroim hoop toision in conical tank. 


( 5 ) 



( 6 ) 

(7) 
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rnaxiniijin deflection of th’e tower/ and 

maximuin ultimate resistance of a local section of unit 
width of supporting shaft at its bottom in case of 300 
kl capacity tank. 

The resistance of the structure to collapse loads is 
found out to be more than adequate/ in all the cases, 

Ihus the serviceability limit state of cracking is the 
most critical one in case of reservoir portion/ while limit 
states of deflection and ultimate strength of a local section 
of the supporting shaft govern its design. 

3.10 Method of Using the Design Tables 

3,10,1 Example 

The method of using the design tables can be illustra- 
ted by the following example: 

A conical water tank supported on a cylindrical ^aft/ 

3 

is reqtiired to be designed for a capacity of 300 m , 

The staging height is 25 m and the wind and seismic 
zones to be considered are II and IV respectively. 

The Table 3,21 is the relevant design table for 300 

3 

m capacity tank. It gives optimal radius of the supporting 
shaft equal to 1.625 m/ optimal value of ^ 1/^2 ~ ^*0 the 
length of the tank wall L = 8.087 m. The other design vari- 
ables are as follws: 

Thickness of the tank wall * t^’ *= 132 mm 

Areas of hoop reinforcemoit in the three regions (Figure 3,15) 
of the cchical porticai (as calculated from the percentages of 



reinforcement given in the design table) are 1894^ 1438 and 
892 

In the conical portion a minimum reinforcement = 330 mm /m, 
in the meridional direction, is found to be adequate for all 
the cases except for a ^ort length of 1 m near junction 2 in 
which the reinforcement (as per the percentage of reinforce- 
ment given by *Prc 4 ‘ ) ^or each case will be different. For 

2 

this case it is = 660 nm /m. 

Thickness of the supporting shaft 't^' = 143 mm 

O 

Areas of vertical reinforcement in the three regions of supp- 
orting ^aft (Figure 3.16) are 636, 1763 and 2699 mm^m. 

Area of circtamferential reinforcement for supporting shaft 

2 

corresponds to the minimum is equal to 357.5 ran /m, 

2 

A 75 ran thick top dome with minimxjm reinforcement of 187.5 ran 

in both the directions will be adequate for all the cases. 

Whereas the inner cylindrical shaft of 125 mm thickness with 

2 

minimum reinforcement ** 312,5 mm /m, in both the directions, 
is found to be adequate for all the cases except for a short 
length of 1 m near junction 2 in vdnidh the vertical reinforce- 
ment for each case will be different. For this case it 
2 

is *= 750 ran /m. 

Ring beams of 200 x 200 mm size with 0.50% hoop reinforcement 
are found to be more than adequate at both the junctions , for 


all the cases 



CHAPTER 4 


OPTIMAL DESIGN OF TYPE II WATER TOWERS 

4.1 Design Considerations 

For 500 to iOOO kl capacity tanks, water towers of 
the type shown in Figure 1.1(b) are preferred. In this type 
of towers, conical tank is provided with a spherical bottom 
and the diameter of the inner cylinder is kept smaller than 
that of the supporting shaft. This arrangement results in 
more space for storage of water as compared to Type I water 
towers and hence leads to overall economy. The diameter of 
the inner cylinder is fixed from the consideration of minimum 
space required for an easy access and in the presant study it 
is kept equal to 1.5 m. 

The size of the opening of the bottom dome is quite 
large as compared to its overall dimension as a result of 
which the influence of edge-loadings at two ends interact and 
both membrane arid banding forces are found throughout the dome 
Although the semi-central angle of this dome can be chosen as 
design variable, its cost in relation to the cost of the 
system is very small. Therefcare, a suitable value which 
miniroisBes the design forces in the dome is obtained indir- 
ectly , It is found by actual analysis that the magnittides 
of banding moments become minimal for 45® (0^^ as 

defined in Figure 3.6), in most of the cases. Hence a value 
of 45® is chosen for 0 ^ for the bottom dome. As the siae of 
the bottom dome is small ani that too with a large opening 
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at the top, the whole dome becomes a critical region wherein 
stress variations are violent. Keeping this in view, as a 
precautionary measure, minimum thickness of the dome is 
restricted to 150 mm with reinforcement not less than 0.5% in 
both the directions, Ihe reinforcement is equally distributed 
on each face in a staggered fashion. In addition to these, 
all the design considerations discussed in Secticn 3,4 will 
also be applicable to the design of these towers. However, 
the thiclcness of the top dome is kept 100 mm, instead of 75 
mm as the design forces will be more because of the larger 
size of the dome as compared to that required for Type I towers . 

4,2 Elastic Analysis 

In order to find out the design forces in the reser- 
voir porticm, elastic analysis is carried out using the force 
method of analysis already described in Section 3.3,1. The 
Figure 4.1 shows the water tower under consideration, vAiereas 
Figure 4.2 shows the same, separated into simple shell elements 
with corrective loadings at various junctions. 

The compatibility of deformations and equilibrium | 

conditions at junction 1 yield an equation which will be | 

exactly the same as giv«i by the Eq, (3,16). Whereas, for 

: 

junction 2, these conditions will result in an equation of | 

the fcanm giv«a by Eq, (4.1). s 

■ I 




zao 
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( 4 . 1 ) 
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where 


oBD2 , r,BD2 

Op and Pp are the displacemait and rotation of 

bottom dome at junction 2 due to primary loadings . 
,BD2 ,BD2 ^BD2 ^ HD 2 

*■11 ' ^12 * ^21 ^22 flexibility coeffi- 

cients which represent displacemoits and rotations of 
bottom dome at junction 2 due to unit edge loadings 
on the dome at that junction itself, and 


cxii , ct^2 ' 


and cf 22 ^ are the coefficients 


which give displacemaits and rotations of bottom dome 
at junction 2 due to unit edge-loadings at the other 
end of the dome. 


Other terms have already been defined in Section 3. 3. 1,2. 


Similarly the compatibility and equilibrium conditions 
for junction 4 will give a set of linear algebraic equations 
which can be expressed by Eq. (4.2). 

Now Eqs, (3,16), (4.1) and (4.2) are combined together 
to form a single set of twenty linear algebraic equations in 
twenty un)cnowns to arrive at the solution. The various flexi- 
bility coefficients and the deformations due to primary load- 
ings, required for the solution, can be obtained from Tables 
3,4 through 3.16 and 4.1, 

The procedure given in Section 3.3.2 is adopted to 
catry out ■tdie elastic analysis for finding the design forces 
in the supporting shaft. 


4.3 Limit Analysis 

The limit analysis of conical tank discussed in 
Section 3 . 5.2 will be valid for this type of towers as well. 



ib:d 



( 4 . 2 ) 




(c) Equally distributed load 'p’ per unit length 
the edge 

Fig. 4-3 Primary loadings on bottom dome 



Table 4.1 Membrane solutions for opoi spberlcal shell - primary loadings 
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Itie possibility of collapse of the bottom dome prior to the 
ccwiical tank is also rare as it is primarily subjected to 
membrane compression* The limit analysis of the supporting 
shaft given in Section 3.5.3 will be applicable here also, 

4.4 Optimization 

4.4.1 Design variables 

In addition to the variables considered for the 
Type I tower * two more design variables / one for the thick- 
ness, and the another for the percentage of reinforce- 
ment, bottom dome are considered in the optimal 

design of Type II water towers. Thus, the total number of 
design variables become twelve for this case, 

4.4.2 ObJLective function 

The objective function of such a tower includes the 
cost of concrete, reinforcement and formwork for the bottom 
dote in additican to the costs considered for finding the 
objective fiinction in case of Type I towers, ihus, the 
objective function can be writtai in the form 


F » 5‘V + 2 it R^(cosj2f2bd " ^°®^^lbd^^^ 


^ ^ ^rbd 


tj^ Cs t 2Cf ) 


(4.3) 


vAiere F‘ is the value of objective function given by Eq, 
(3,39) for a tower of Type I, ^ 

a radius of the bottom dome. 
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= thickness of the bottom dome^ 

^rbd ~ peJ^centage of reinforcement in the bottom dome 
for both circumferential and meridional direc- 
tion # and 

^2bd ^Ibd corresponding values of 0^ 

(Figure 3,6)/ for the bottom dome, 

4.4.3 Constraint s 


The first thirty constraints for this type of tower 

* 

will be essentially same as those considered for the Type I 
towers. In addition to these, the constraints governing the 
design of the bottom dome are needed to be incorporated. 
Serviceability limit state of cracking demands that the flex- 
ural tensile stresses in concrete as well as in steel should 
remain within prescribed values as stated in Section 3,6.4. 
Ihus, if and are the maximum bending moment carrying 

capacities of unit width section of the bottom dome as per the 
limiting stresses for concrete and steel respectively, then 
the constraints thirty-one and thirty- two can be expressed in 
the form 


s 

|s 

i; 

j; 


- 1 < 0 , (4.4) 

”bdc 


and S ° <■‘•5) 

“bds 


where 


* 


actual nraximum balding moment 


(maximum of 


M,y and M^) per ijnit width of the bottom dome, 
at service loads. 
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Similarly for safety, the limit state of strength 

requires that the ultimate mcsnent of resistance 'M , per 

urbd ^ 

width of bottcxB dome should be greater than or equal to 

the actual value of maximum ultimate moment ‘M , . Thus, 

uabd ' 

the thiirty— third constraint can be writtoi as 


uabd 

M 

urbd 


1 < 0 


(4.6) 


In addition to these, two more constraints, regarding 
the minimum thickness and reinforcement for the dome can be 
ejqaressed as 


^min 

%d 


< 0 


(4.7) 


min 

Pbd 


1 < 0 


(4.8) 


Thus, the total number of constraints, for this case, ccames 
out to be thirty -five. 

4.5 Parametric study 

For water towers of Type II, optimal limit state 

design has been carried out for the following cases: 

3 

(1) capacity in ra : 500, 750 and 1000; 

(2) staging heis^t in ms 15, 20 and 25; and 

2 

(3) basic wind pressure in WN/n : 1.0, 1.5 and 2.0. 

A free board of about 200 ram. has been provided in all 
tiie cases. The sxoitability of these designs also, obtained 
lor differ^t wind zones, is established for the various 
seismic zcmes of India. 



Table 4.2 Optimal designs of Type II towers with 500 m capacity tank 
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'able 4.6 Optimal designs of Type II towers with 750 m capacity tank for non-optimal shaft 
radii 
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Table 4,6 (continued) 
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l^ble 4,7 Optimal designs of Type II towers with 1000 m capacity tank for non-optimal shaft 
radii 
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Table 4.7 (continued) 
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Table 4*8 


Effect of cost ratios on optimal configuration 
Capacity of tank = 750 L = 25 m, pj^ = 1.0 kN/hi' 


s. 

No. 

CRl 

CR2 

3^ 7i 

w to 

R/Rj 

Optimal 

cost 

(Rs) 

Optimal 

configuration 

1 
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332490 
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= 1.75 m. 
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X z 
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a = 1.020 rad 
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0.10 
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R /R = 4*5 or 
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0.15 

1.5 

4.5 
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5.5 
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R^ = 1.75 m. 




1.75 
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R,/R^ * 4,5 or 





5.0 

373066 

X z 




2.0 
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370837 

a = 1,020 rad 





4.0 

362330 






4.5 

384219 
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Table 4.9 Effect of cost ratios on optional design variables 
Capacity of tank = 750 L = 25 m, Pj^ = 1-0 m/m' 


R, {CRl K:R2 

-*« 1 h 


Optimal design variables 


^c 

^rcl 

r ""'""" 

^rc2 

"'""1 

^rc3 

r -* : 

^rc4 

Prbd 

(nsm) 

i 

1 

^rsl 

^rs2 

Prs3 

P 

^rcy 

^bd 

(ram) 


Optima 

cost 


1 1.75 4.5 1.25 0,10 223 1.459 1.060 0.663 0.25 0.50 314087 

214 0.258 0.585 1,324 0.25 150 


2 1.75 4.5 1.25 0.15 219 1.450 1.083 0.708 0.25 0.50 357754| 

210 0.292 0.674 1.236 0.25 150 i 


3 1.75 4.5 2,00 0.10 222 1.405 1.070 0.670 0.25 0.50 390498 

217 0.280 0.658 1.158 0.25 150 

4 1.75 4.5 0.75 0.10 219 1.454 1.085 0.681 0.25 0.50 260065 


213 0.262 0.742 1.115 0.25 150 


1 
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Table 4.10 Comparative optimal cost study of Type I and 
Type II water towers 

L « 20 m, p^ = 2.0 kN/m^ 


Tank 

capacity 

(m^) 

Optimal cost of 

Type I tower I 

1 

(Rs) i 

I 

Optimal cost of 
Type II tower’ 

(Rs) 

I Difference in 
i costs of two 

1 towers 

1 (Rs) 

100 

99196 

97197 

1999 

200 

132777 

130799 

1978 

300 

164292 

158762 

5530 

500 

245902 

209878 

36024 

750 

347775 

296154 

51621 

1000 

456750 

379031 

77719 
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4,6 Results and Discussions 

Similar to the case of Type I towers / for this type 
of towers also it has been observed that for optimal configu- 
ration the bending moment and shearing force in the reservoir 
portion reduce considerably and the value of maximum hoop 
tension becomes the minimum/ in most of the cases, 

'Ihe Tables 4,2 through 4,4 give optimal designs of 
Type II towers for various capacities/ staging heights and 
wind and seismic zones. Whereas, Tables 4.5 through 4,7 
give optimal designs for these cases when the radius of the 
supporting shaft is other than the optimal one but with a 
value in its neighbourhood. 

Itom the results given in Tables 4,8 and 4,9 it can 
be observed that, for this type of towers also, the effect of 
variotis cost ratios on optimal configuration and design vari- 
ables is quite similar to that already found for Type I towers. 
Thus, optimal ccxifiguration remains iinaffected by any change 
in the cost of materials. However, a slight variation obser- 
ved in the optimal values of some of the design variables is 
negligible for all practical purposes. Hence pres«it optimal 
designs will remain optimal or almost optimal, in future, 
for different costs of materials. 

The optimal value of a is found to lie between 0,991 
to 1,111 radian for all ihe cases but for roost of these it 
is between 0,991 to 1.02 radian. The optimal value of is 
found to be 1,75 m for most capacities, staging heights and 
wind and seismic zones. 



Prom the Table 4,10 it is evident that for capacities 
upto 300 kl the difference between the costs of two types of 
towers viz. Type I and Type II is marginal (1 to 3.5%). How- 
ever for capacities 500 kl and above the Type II towers are 
quite economical » Although the cost of the Type I towers for 
capacities upto 300 kl is slightly higher than the Type II 
towers^ from the constructional point of view the former 
are simpler , Hence the preference of Type I towers over 
Type II upto 300 kl capacity is justifiable. 

Discussion on relative values of seismic and wind 
forces as well as those regarding constraints which are active 
except the ones dealing with minimum thickness of tank wall 
and supporting shaft given in Section 3.9, is found to be 
valid for this type of towers also. In addition to these, 
the constraint which deals with the ultimate strength of a 
section of the conical shell at junction 2 is also found to 

3 

be active in some cases of 1000 m capacity tank. 

4.7 Variable Ihickness of the Conical Tank Wall 

Althoi:^fh the analysis as well as the optimal designs 
of the water towers described in the foregoing are strictly 
applicable when the thicknesses of the various elements are 
unifo»m 4 f still the thickness of the wall of the conical tank 
can be varied without causing any serious error in the 
values of the forces ^ch govern the design of the tank. 

The location of the maximum hoop tension in the conical tank 
is invariably found to be betweeii 0.3 L to 0,5 l from the 
junction 2, for all the cases. From 0,5 L to the junction 1, 



the hoop tension decreages uniformly and becomes almost equal 
to zero at Junction 1, Thus, the thickness of the ocsaipal 
tank wall can be decreased tiniformly from the maximum, equal 
to that obtained in the optimal design, at 0*5 1 to a minimum 
prescribed thickness, 125—150 mm, at the Junction 1. The area 
of hoop reinforcement to be provided in each of the three 
regions will remain equal to that provided for the optimal 
design. This variation of the thickness will not lead to any 
error , particularly , in the values of hoop tension at loca*^ 
tions which govern the design of the tank wall and also does 
not violate any of the governing design ccaistraints. However, 
there may be a little change in the bending moments in the 
edge region of Junction 1 which can be taken care by providing 
slightly thicker section for a short length near junction 1. 
Such a taper results in saving of a considerable volxime of 

3 

concrete, particularly, in case of 750 and 1000 m capacity 
tanks. 

4,8 Example 

To illustrate the method of using the design tables 
for the type of the towers considered, a design example has 
been tak«i. 

Problems A conical water tank supported on cylindrical 

shaft is required to be designed for a capacity of 
750 staging hei^t 20 m and wind and seismic 
zones II and III respectively. 

' 3 

The Table 4.3 gives the optimal designs for 750 m 
capacity tanks , From this table the various design variables 



184 


can be found for the given staging height and wind and seismic 
zones. The optimal design values are as follows: 

Radius of the supporting shaft 'R^' = 1.75 m. 

Value Of ~ 4.5 (a = 1.020). 

Length of the ccnical wall *L* = 11.70 m. 

'Ihictoess of conical tank wall ‘t ' = 221 mm. 

c 

Areas of hoop reinforcement in the three regions of the conical 
portion (Figure 3.15) in order are 3158^ 2358 and 1480 mm^/m. 
Area of meriodimal reinforcement in the conical portion (as 
per minimum requirement of 0.25%) = 552.5 mmVtn* 

Ihickness of the supporting shaft ' t_' = 231 nm. 

Areas of vertical reinforcement in the three regions of the 
supporting shaft (Figure 3,16) in order are 649, 1257 and 
2136 mmVm. 

Area of circiunferential reinforcement for supporting shaft (as 

2 

per minimum requirement of 0.25%) = 578 ram /m. 

Hiiokness Of the bottom dome • t^' = 150 mm 

Area of reinforcement in meridional as well as in circumferen- 

2 

tial direction for the bottom dome is = 750 nan /in. 

2 

A 100 mm thick top dome with 250 mm /m reinforcement in both 
the direction*and the inner cylinder of 125 mm thickness with 
312.5 mmVm reinforcement in the two directions is found to 

be adequate# for all the cases. 

Ring beams of 330 x 330 mm and 225 x 225 rran sizes 
with 0*5% hoop reinforcements will be adequate enough for 
junctims 1 and 2# and 4 respectively. 
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Ihe variation of the thickness of the wall of the 
conical tank can be done as per the procedure described in 
Section 4.7. The thickness of the wall from junction 2 to 
0.5 L (5.85 m) will remain equal to 221 mm itself. From 0.5 L 
(5.85 m) to the point of junction 1 it can be uniformly decre- 
ased to a minimum of 125-150 mm. The area of hoop reinforce- 
ment will essentially remain unchanged in all the three regions 
and will be equal to 3158,. 2358 and 1480 mm^m respectively. 
However, for a short length near junction 1 a slightly thicker ' 
section with gradually increasing thickness should be provided. 



CHAPTER 5 


OPTIMAL DESIGN OF TYPE III WATER TOWERS 


5.1 Design Considerations 

Scxnetimes a reservoir is reqxiired to be provided 
with a partition for various reasons. In this chapter, water 
towers in which reservoirs having partitions are considered. 

Hie Figure 1.1(c) i^ows a tower of this type having a cylin- 
drical partition which divides the reservoir in two portions. 
Water towers of this type are also considered for 500 to lOOO kl 
capacities. 

Design forces for the reservoir portion are obtained 
by ccaisidering all the possible governing axisyrametric loading 
combinations such as, both the portions of the tank are full, 
only outer portion is full and when only inner portion is 
full. It has been observed that the bending forces in edge 
region of junction 2, in all the elements, are found to be 
more when only outer portion is tvll, whereas the governing 
design forces, particxilarly the hoop tension, for the cylin- 
drical partition is obtained when only inner portion is 
full (Figure 5.3). In addition to these, all the design 
considerations which are discussed in the Sections 3,4 and 

4.1 will be valid for this type of towers also. The minima 
thickness of the partition wall is restricted to 12S imn. 
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5.2 Elastic Analysis 

The force method of analysis described in Section 
3.3,1 is used for elastic analysis for various axisymmetric 
loadir^ combinations discussed in the previous section. The 
Figure 5,1 shows the details of the water tower and the 
Figxare 5.2 shows the same, exploded into simple shell elements 
with corrective edge— loadings at different junctions. 

At junctions 1 and 4, the compatibility of deforma- 
tions and equilibrixim conditions give relations which will be 
exactly the same as those given by Eqs, (3.16) and (4.2). 
Whereas, for junction 2 it will be slightly different from 
the one givan for Type II towers. For this type of towers, 
equilibrium and compatibility conditions at junction 2 will 
yield an equation of the form given by Eq. (5.1). in which 

ill , f ^2 / ^21 ^22 appropriate flexi- 

bility coefficients for the cylindrical partition, 
and 

represent displacanent and rotation 
respectively, of the cylindrical partiticxi due to 
primary loading, at juncticm 2. 

The flexibility coefficients and deformations due to primary 
loading for cylindrical partition can be obtained from Tables 
3.6 and 3.10. The results of the Table 3.10 should be vised 
with proper sign. These results are to be used with same 
sign when hydrostatic pressure causes hoop compression and 
with oppcxiite sign in case of hoop tension. 
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Thus the Eqs, (3,16)^ (4.2) and (5.1)# when combined 
together# form a single set of twenty-two algebraic eqxiations 
in twenty-two unknowns (X1-X22) which can be solved by matrix 
inversion. The various elements of the flexibility matrix 
and the components of the deformation vector due to primary 
loading can be obtained by using the results given in Tables 
3.4 through 3.16 and 4,1,. 

For finding the design forces in the supporting shaft#^ 
the procedure given in Secticxi 3,3.2 is adopted# for the elas- 
tic analysis. 

5.3 Distribution and Arrangement of Reinforcements in the 
Cylindrical Partition 

The Figure 5.3 shows general pattern of variation of 
stress resultants along x-direction for the partition. It 
can be seen that the value of hoop tension 'Nq' is small for 
a very short length at the bottcxn of the partition and it 
shoots up to its maximum value at ^ 0»1 (where 1 is the 

length of the cylindrical partition) and then gradually decre- 
ases to zero value at its top. Cxortailment of hoop reinfor- 
cement for such a short length at the bottcxn will not be of 
imich use# therefore# a stepped distribution of reinforcement 
with points of curtailment at one-third and two-third heights 
of the partition has been adopted in all the cases. Simi- 
larly# the maximum bending moment occurs at the bottom 

and it decreases at exponential rate. Therefore# a minimu^^^ 
reinforcement of 0*25% will be adequate in x-direction except 
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than 1 m in all the cases considered in the present study. 
Hence in finding the objective function/ this length is taken 
as 1 m. The hoop reinforcement as well as the reinforcement 
in the vertical direction is provided on both the faces, 50% 
on each side, in a staggered fashion. 

5.4 Limit Analysis 

The limit analysis of conical tank and that of supp- 
oirting shaft described in Sections 3,5.2 and 3.5.3 will be 
valid for this type of towers as well. In addition to these, 
there is a possibility of collapse of cylindrical partition 
also when only inner portion of the tank is full, in the 
limit analysis of structures, it is a common practice to get 
a lower bound on the load carrying capacity from a statically 
admissible stress field or to obtain an upper bound from a 
kin«natically admissible mechanian. In case of reinforced 
concrete shells it is very much tedious to obtain an exact 
solution through lower bound approach except for some simple 
cases. However, in the case of cylindrical partition this 
can be applied without much of difficxxLty, For this, the 
solution commences from the eq\iilibri\im equation and a kine- 
matically admissible mechanism for the stress field satis- 
fydng the equati<xv of equilibrixim is then found. Hence the 
solution so obtained is complete and gives the exact limit 
load. Sawczuk and Olszak (1961) have obtained the criterion 
of failure for an axisymraetrically loaded cylindrical shell 
of reinforced concrete, which is briefly described in the 
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5.4.1 Criterion of failure 

For an axisymmetric shell there are, in gaieral, four 
stress resultants; the bending moment M , the shear force Q . 
and membrane forces N and N , When Q is neglected in the 
condition of failure, the failure criterion is represented by 
two parabolic cylinders bounded by two parallel planes. If 
force N is small in comparison with N_ and M , which is tme 

•A. y X 

for the case under consideration in the present study, the 
intersection of the cylinders with the plane N =0, gives the 
simplified failxire criterion in the M -N plane. The Figure 
5.4 shows the simplified failure criterion and the rule of 
instantaneous deformations. In this figure, and X denote 
the generalised strains associated with and Ng and = 
ultimate moment of resistance per xmit width, to resist posi- 
tive moment; M' = ultimate moment of resistance per unit width, 
o 

to resist negative moment; = ultimate hoop compression 
capacity per unit width; and *= ultimate hoop tension capa- 
city per unit width, all pertaining to ariy section of the 
shell under consideration. The moment and membrane force 
capacities can vary from, section to section, but the profile 
of the failure criterion remains essentially the same. 

The Figure 5.5 shows the forces acting on an element 
of the shell under consideration and the sign convention 
adopted. The equation of equilibrium for this case of 
cylindrical shell can be written as 

.= 0 ( 5 . 2 ) 
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where = intensity of hydrostatic pressure, 

R = radi\xs of the cylindrical shell, and 
w = radial displacanent. 

The generalised strain rates can be expressed in 

terms of w in the form - d^^/dx^ and = v^/r. From the 

Pigvire 5.4 it can be seen that regimes AB and CD have to be 

rejected because along AB and CD, ^ is equal to zero and 

2 . 2 

d w/dx is not equal to zero, which is impossible. Along BC 
and AD, is equal to zero which leads to possible collapse 
mechanism of the form w = c^x + C 2 » Therefore the collapse 
mode is always conical in form. The exact shape of 
the collapse mechanism depends on the shell parameters, and 
moment and hoop tension capacity distribution provided. 

Based on these, shells may be classified as short, medium and 
long, 

5.4.2 Collapse mechanism 

Ihe size of the cylindrical partition, considered in 
the curr«iit investigation is such that it will be always in 
the category of long cylinders for which the failxire will be 
caused by partial collapse of the bottom portion of the shell, 
in all the cases. Details of this mode of collapse are 
given in the Figure 5,6. Positive hinge circles form at x = 0 
and X * X 2 and a negative hinge circle at x = x^. Beyond 
X * X 2 # the bCTKiing mom^t diagram can be continued in a 
statically admissible manner. Ihe ultimate hoop tension 
• . .. -1 4 1 on conforms to the hoop reinforcement 
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distxibution ^ described in Section 5 . 3 / provided for crack 
control at service loads. The value of ultimate hoop tension 
capacity jiamps from to N ^2 the second region bovindary 
X = I 2 and from N ^2 to at the first region boundary x = 1 ^/ 
which is permissible; but the bending moment and ^ear force 
have to be continuous. Due to the hydrostatic loading the 
value of p^ will be equal to p^Cl - ^) / where p^ is the 
collapse pressure. The equation of equilibrium can be written 
as ■ ■ ■ ' 


d^ 


N. 


dx' 


X 0^/1 Xv 

7 ^ - R Pq(1 - J-) 


(5.3) 


Integra tioi of this equation, in the bottom two 
regions, results in four .constants of integration, in addi- 
tion to the moment and shear continviity cc»iditions at the 

first region boundary x = 1 , , the condition that M = M-., at 
dM ^ X u± 

X * 0 and = 0 at x = X 2 are used to evaluate these cons- 
tants of Integration, The value of these constants will 
depend upon the location of third hinge circle. For the type 
of hoop tension capacity distribution considered herein, the 
third hinge will usually be in the bottom region or in some 
cases it may be in the middle region. Depending upon the 
locaticxi of third and second hinge circles, three cases are 
considered as under ! 

Case I 


When the third hinge circle is assumed to be in the 
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simplification^ can be given by the following equations: 


dM p 

X „ -^o 


N. 


and M. 


(21x - x^ - 21 x 2 + X 2 ) - -|i(x - X 2 ) , (5.4) 

P. 


X 


gY (31x^ _ ^ 


6IX2X + 3X2X) 


N 


tl / 2 

(x - 2 X 2 X) + Mq^ 


2R 


(5.5) 


In order to solve for the three unknowns x, and 

1 

X 2 / the following conditions are used: 


at 

X 


(region 1) 

11 

0 

(5.6) 

at 

X 

=* X^ , 

(region 1) 

“x = -%2 ’ 

(5.7) 

at 

X 

= X 2 , 

(regicxi 1) 

“x = “03 

(5.8) 


Case II 

In some cases, the third hinge circle may form in the 
middle region whereas the second one remains in the bottom 
region itself, For this case expressions for shear and 
moment are required for both bottom and middle regions. The 
simplified equations for Shear force and bending moment in 
the two regions are the following: 
dM 


region Is *= - x - 21x2 + X2) 

-(1, - Xo) 


\2 

-^(x-1/)- ^ 

R 1 


R '*1 ”■ ■^2'' ' 


(5.9) 
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M 


X 


region 2 j 


and 


dM 

• _ > 

dx 


M 


P© 2 3 2 

^(31x - X - 6 IX 2 X 4- Sx^x) 


2 ^t2 

2r (x - 21j^x) - ^ (l^x - X2X) + , 


(5.10) 


2 2 ^t2 

~(21x - X - 21 x 2 + xp - 3 C 2 ) / 


(5.11) 


X 


■g^(31x^ - x^ - 6IX2X + Sx^x) 


N 


t2/ 2 


2R 


(x - 2X2^ ^ 1 ^ 


+ + 


\l^l 


01 2R 


(5.12) 


In order to solve for the three unknowns x^^ and 
X 2 / the following conditions are used: 

dM 



at 

X 

= x^ , 

(region 

1) 

II 

0 

(5.13) 


at 

X 

= X^ , 

(region 

1) 

^x = - %2 ' 

(5.14) 

and 

at 

X 

= X 2 / 

(region 

2) 

^x = ^03 

(5.15) 

Case 

III 



• 





When both second and third hinge circles form in the 
middle region (Nq = - N^ 2 ^' equations for shear force and 
bending moment for the two regions will rapnain exactly the 
same as those given for Case II. using these equations, the 
three unknowns p©, x^ and X 2 can be obtained using the foll- 
owing conditions: 


dM 


at X = Xj^ # (region 2) = 0 , 


(5.16) 
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at X = , (region 2) , (5.17,) 

and at x = X 2 / (region 2) (5.18) 

Application of such conditions, in each of the foregoing cases y 
will result in three nonlinear simultaneous equations which 
have been solved numerically rasing the Newton- Raphson method, 

5.5 Optimization 

5.5.1 Design variables 

k 

All the design variables considered for the Type II 
towers will be valid for this case also. In addition to 
these five more design variables, first one for the thickness, 
tp, of the partition wall, next three for the percentages of 
hoop reinforcement, p^^pj^^r Pj.p 2 and last one for 

the percQitage of reinforcement, Pj^p^/ in the vertical direc- 
tion in the cylindrical partition for a ^ort length equal 
to 1 m near jranction 2. Thus, the total number of design 
variables, to be considered in the optimal design of Type III 
water towers, is seven te^, 

5.5.2 Objective function 

The objective function for the case \inder considera- 
tion will include the cost of materials and formwork required 
for the partition in addition to the costs considered for 
finding the objective function of the Type II towers . Hence 
the expression for the objective function in this case will 
be of the following form: 



F = 


F" + 2 IT R 2 I [t„C^ + 


0.25 


Pe, C„ + 2C, 


p''c ' 100 '■p ''s ''s ' “^f 

'rpl ^rp2 ^rp3 


t p C 

^200 Pt-i-iO Pi-n '^^3 


- 0.25 

+ 2’‘«2 Ps = S > 


(5.19) 


where F" is the value of objective function for the Type II 
towers as given by Eq. (4.3)/ and 1 is the length of the 
cylindrical partition. 


5.5.3 Constraints 

All the thirty-five constraints which are considered 
for the Type II towers will be valid for this case as well. 

In addition to these, the constraints governing the design of 
the cylindrical partition are to be incorporated in the presoit 
case. The serviceability limit state of cracking demands 
that the direct and flexural tensile stresses in concrete as 
well as in steel should remain within prescribed values as 
stated in Section 3.6,4. Thus, if ^©p2 ^©p3 

sant- the maximum hoop tension in regions 1, 2 and 3 respec- 
tively and M and M are the maximum bending moment carrying 
capacities of unit width section of the partition as per the 
limiting stresses in concrete and steel respectively, then 
the constraints thirty-six through forty-three can be written 
as: 


N /ft (-££■-) f I 

©pr ^^P^IOO ' ^st max;;} 


N /ft (-£H.?.) f 
^0p2/*-^p4oO ^ ^st max 


] - 1 < 0 , 

] - 1 j< 0 , 


(5.20) 

(5.21) 
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^eps/^^^p^ioo^^ ^st 

max ^ 

0 

VI 

1 

(5.22) 

M 

P^ i ^ 

pc 

0 , 



(5.23) 

M 

-ES. - 1 < 

Ss - 

0 , 



(5.24) 

NQpi/(Ae^l 

f a. ) 

ct max 

- 1 

< 0 , 

45.25) 

N0p2/(Ag^2 

f ) 

ct max 

- 1 

< 0 , 

(5.26) 

^ep3/^^eu3 

f a. ) 

ct max 

- 1 

< 0 

(5.27.) 


where = actual maximum beriding moment per unit width of 

partition wall in the edge region of junction 2 , 
and 


^st max' *ct max' defined in the 

Section 3.6.4 except in finding t^ and 

Pj.^^ are to be replaced by t^ and Pj.pj^ for the 
present case. 

The constraints regarding the minimum thickness and 
reinforcement for the partition wall can be expressed in the 
following forms 


t . 
min 

P 

mxn 

Prpl 


0 , 



^r min 
^rp2 


1 




C5.28) 

( 5. 29 ) 


< 0 


(5. 30) 





and 


min 

^rp3 


1 < O , 


. 1 < 0 

Prp4 


(5.31) 


(5.32) 


Ihus# for this case^ the total nximber of constraints 


is forty-eight. 


5.6 Parametric Study 


Similar to Type II tower, in this case also, optimal 
limit state design has been carried out for the following cases i 


( 1 ) 

( 2 ) 

(3) 


3 

capacity in m : 500, 750 and 1000; 

staging height in ms 15, 20 and 25; and 

2 

basic wind pressure in Idl./m ; 1.0, 1.5 and 2.0. 


Suitability of these designs are also established for 
various seismic zones of India and a free board of about 200 
mm has been provided in all the cases. 


5.7 Results and Discussions 

Similar to the previous two cases, for this type of 
towers also it is found that for optimal configuration the 
bending moment and shear force reduce considerably in all the 
elements in the edge region of junction 2 and the value of 
maximxom hoop tension in conical tank, becomes minimum, in 
most of the cases. 

The optimal designs of this type of towers can be 
obtained from Tables 5,1 through 5.3, for various capacities, 
staging heists and wind and seismic zones. Whereas, optimal 



Table 5.1 Optimal designs of Type III towers with 500 m capacity tank 
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Table 5,5 (continued) 
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Table 5.6 Optimal designs of Type III towers with 1000 m capacity tank for non-optimal shaft rada 
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Table 

5.7 

Effect 

of cost 

ratios on 

optimal 

configuration 



Capacity of tank = 1000 

L = 

25 m, pj^ = 2.0 WSI/m^ 

S. 

CRl 

CR2 

^2 

R./R^ 

Optimal 

Optimal 

No, 





cost 

configuration 




(m) 


(Rs) 


1 

1.25 

Q.IO 

1.75 

4,5 

500601 






5.0 

474839 






5.5 

477288 

R« = 2,0 




2,0 

4.0 

486777 

z 





4.5 

462523 

R,/R^ = 4.5 or 





5.0 

476437 

X z 




2.25 

3.5 

502807 

a = 0.961 





4.0 

475360 






4.5 

486921 



2 

2.0 

0.10 

1.75 

4.5 

622368 






5.0 

590239 






5.5 

598344 

R^ = 2.0 m. 




2.0 

4.0 

602469 

z 





4.5 

577346 

R 1 /R 2 = 4.5 or 





5.0 

598896 





2,25 

3.5 

621464 

a = 0.961 





4.0 

589094 






4.5 

615327 



3 

1.25 

0.15 

1,75 

4.5 

565562 






5.0 

542790 






5.5 

540032 

R^ = 2.0 m. 




2,0 

4.0 

556466 

z 





4.5 

530409 

Rt/R^ > 4.5 or 





5.0 

546127 

X z 




2.25 

3.5 

566368 

a = 0.961 





4,0 

544627 






4.5 

560431 



Table 5,8 Effect of cost ratios on optimal design variables 
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Table 5.9 Capacity of inner portion of the tank for different 
cases 


Tank 

capacity 

(mP) 

nr 

1 

1 

! (ra) 

1 

* - 

T — r 

1 1 

; a ; 

1 • 

1 1 

> 1 

J (radian) ; 

• 1 

1 ..-- - 1 

1 

Height of ; 

partition wall,* 

1 1 

(m) j 

1 

Capacity of the 
inner portion 

(ra^) 

500 

1.50 

0.908 

7.679 

40.71 


1.75 

0.837 

6.792 

53.34 



1.013 

8.410 

66.05 


2.00 

0.817 

6.397 

69.08 



1.018 

8.109 

87.57 

750 

1, 50 

0.960 

9.641 

51.11 



1.090 

11.510 

61.03 


1.75 

0.875 

8.388 

65.88 



1.029 

10,173 

79.90 


2.00 

1.011 

9.571 

103.36 


2.25 

0.840 

7.530 

106.45 


1.50 

1.075 

12.490 

66.21 

1.75 

0.998 

10.864 

85.32 

2.00 

0.961 

10.020 

108.21 

2.25 

0.965 

9.740 

137.70 


1000 
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designs for these cases for non-optimal radii of the supporting 
shaft are given in Tables 5.4 through 5.6. It can be observed 
fron the resiilts that the optimal value of and a for a 
particular capacity, staging height, and wind and seismic 
zones, for this case, is found to be less as compared to that 
of Type II towers. This is due to the fact that the cost of 
cylindrical partition reduces with decrease in the value of R 2 
and a both. 

For this type of towers the optimal value of a lies 

between 0-817 to 1.090 radian in all the cases but for most of 

these it is between 0.9 to 1.03 radian. The optimal value of 

3 

R 2 for most of the cases of 500 and 750 m capacity tanks is 

3 

1.5 m, whereas it is 1.75 m for the towers with 1000 m capa- 
city tank. 

Discussion on relative values of wind and seismic 
forces for various cases, effect of cost ratios on optimal 
configuration and the design variables as well as those dealing 
with constraints which are active, given in Sections 3-9 and 
4.6, is found to be valid for this type of towers also. In 
addition to these, the following constraints, all pertaining 
to partition wall, are also critically satisfied in the 
optimal design of this type of towers; 

(1) the constraint dealing with minimum thickness, 

(2) all the constraints corresponding to maximum allowable 
stress in hoop reinforcement, 

the constraint which deals with maximum allowable direct 
tensile stress in ccoicrete in the region of maximum hoop 
tensicai, and 


( 3 ) 
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(4) the constraint in relation to maxiraim allowable flexural 
tensile stress in steel in the boundary region of 
junction 2, 

Ihe capacity of the inner portion of the tank depends 

upon the value of R 2 and a. The Table 5,9 gives this for 

different capacity/ R 2 and a values. It can be seen that the 

% 

capacity of inner portion for 500 / 750 and 1000 m capacity 

tanks lies between 40.71 to 87.57/ 51.11 to 106,45 and 66,21 
3 

to 137.70 in respectively for different values of and a in 
each case. For a design case in which there is no specific 
requirement regarding the capacity of inner portion/ the rele- 
vant optimal designs/ from one of the Tables 5,1 through 5,3/ 
can be picked up. When the design of a water tower with a 
specific capacity of the inner portion/ which lies between 
the range given above/ is to be carried out/ then a suitable 
design with optimal value of R 2 or otherwise can be obtained 
from Tables 5,1 through 5,6. However/ these designs are not 
applicable to the cases in which the required capacity of the 
inner portion does not lie between the corresponding range 
prescribed in the foregoing discussion. 

5.8 Variable Thickness of the Conical Tank Wall 

Similar to the case of Type II towers thickness of 
the conical tank wall can be varied without any error in the 
governing design forces and without violating the design const- 
raints. The discussion and procedure regarding variation of 
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thickness of the tank wall/ described in Section 4,7// will 
also be valid for this type of towers, 

5,9 Example 

Similar to the other two cases / to illustrate the use 
of design tables, an example has been considered for this 
case also, 

Problero: A conical water tank supported on cylindrical 

shaft is required to be designed for a capacity 

3 

of 1000 m , staging hei^t 25 m and wind and 
seismic zones I and II respectively, when the 
capacity of the inner portion of the tank is 
required to be not less than (a) 75 m and 
(b) 100 m^, 

3 

(a) It can be seen frcrni Table 5,9 that for 1000 m 

3 

capacity tank, the required capacity of 75 m , of the inner 
portion is available for three values of which are 1.75, 
2,0 and 2,25 m. However, out of these three, R 2 = 1,75 m 
corresponds to the optimal radius of the supporting shaft 
for the giv«i capacity, staging height and the wind and 
seismic zones (Table 5,3), Hence the optimal design values 
obtained from the Table 5,3 are as follows: 

Radius of the supporting shaft '^ 2 ' ~ 1,75m, 

Value of R^/R 2 = 5.0 (a = 0,998) , 

Length of the conical tank wall 'L' = 12,92 m. 

Thickness of the conical tank wall 't^' = 268 mm. 
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Areas of hoop reinforcement in the three regions of the 
conical portion (Figure 3.15) in order are 3722, 2881 and 
1953 

Area of meridicxial reinforcement in the conical portion (as 

2 

per minimvim requirement of 0,25%) = 670 mm /m. 

Thickness of the supporting shaft 't ‘ = 289 mm. 

Areas of vertical reinforcement in the three regions of the 
supporting shaft (Figure 3.16) in order are 803, 1346 and 
2725 mraVm. 

Area of circumferential reinforcement in the supporting ^aft 

2 

(as per the minimiim requirement of 0,25%) = 722.5 mm /m. 

Thickness of the bottom dome = 200 mm.. 

Area of reinforcement in meriodional as well as in circum- 

2 

ferential directicxi, for the bottom dome, is = 1000 mm /m. 

Thickness of the cylindrical partition wall ‘ t^' = 130 mm. 

Areas of hoop reinforcement in the three regions of the parti- 

' 2 j 

tion wall (Section 5.3) in order are 1755, 1397 and 698 mm /m. 
Area of vertical reinforcement in the partition wall for a 
short length of 1 m (as per perrcentage of reinforcement 
•p^p^') = 2639 mm%. 

Area of vertical reinforcement in the partition wall for the 1 

remaining length (as per minimtjm requirement of 0.25%) = 

2 

325 mm /m, 

A 100 mm thick top dome with minimum reinforcement 
of 250 mmVro in both the directions will be adequate for all 
the cases. The inner cylindrical shaft of 125 ram thickness | 
with minimvirm reinforcement = 312,5 mm /m, in the two directions 
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is found to be adequate for all the cases except for a very 
short length of 1 m near junction 2 in which the vertical 
reinforcemoit (as per the percentage of reinforcement given 
) for each case will be different. For this case 
it is = 1375 

Ring beams of 400 x 400 mm and 225 x 225 ram sizes 
with 0-,5% hoop reinforc«nent will be adequate enouch for 
junctions 1 and 2 , and 4 respectively.. 

Ihe thickness of the conical tank wall can be var- 
ied as per the procedure described in Section 4,7. ihe 
thickness from jxinction 2 to 0.5 L (6..46 m) will ranain eqiaal 
to 268 mm itself., Fran this point to the point of junction 1/ 
the thickness can be decreased .uniformly to a minimum of 
125-150 mm, Ihe area of hoop reinforcement in each of three 

regions of the conical portion will essentially remain unchan- 

2 

ged and will be eqiaal to 3722, 2881 and 1953 mm /m respec- 
tively, However, a slightly thicker section with gradually 
increasing thickness for a short length near junction 1 
should be provided, 

(b) In this case the capacity of inner portion of 

3 

the tank is required to be at least 100 m , for which possible 
radii of the supporting shaft are 2,0 and 2,25 m (Table 5.9)., 
However, the optimal design with R 2 = 2,25 ra is not givoi 
for this particular case as it will be very much uneconaiu.— 
cal. Hence the only value of R 2 which will provide the 
required capacity of inner portion is = 2,0 m. This value 
of R« does not correspond to the optimal radius for the case 
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xjnder consideration, Iherefore/ the optimal designs for 
non-optimal supporting shaft radii, given in Table 5.6, will 
be adopted for this case. Adopting the foregoing procedure, 
the design details of the tower can be obtained without any 
difficulty for this case also. 



CHAPIER 6 


COMPUTATICXJAL ASPECTS 

6.1 General 

The developments of modern science have ccaifronted 
the scientist and the engineer with a variety of problems 
which are best solved by approximate methods. Mathematical 
programming techniques/ with the exception of classical 
methods, involve lot of numerical computations. Numerical 
analysis has been applied to scientific and technological 
problems from the very beginning of applied science, but has 
been given new impetus by the development of electronic digi- 
tal computer. In this chapter, ccxnputational details of the 
optimization techniques used in this study, organisation of 
the computer programs, and related aspects are preserxted. 

All the computations have been carried out on DEC-1090 system 
at the Indian Institute of Technology , Kanpur . 

6.2 Algorithms for the Optimization Techniques 

The formulated constrained nonlinear programming 

problems of the present study have been transformed into 

% 

eqxiivalent unconstrained nonlinear programming problems 
through the use of interior penalty fxanctions as described 
in Section 1,4.5. In order to find the minimum of the 
modified objective functicn, Davidcn-FletCher-Powell method 
(DEF method) is used to find search direction S and cubic 
interpolation method is employed to find the appropriate 
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step laigth a« The various steps involved in the solution 
of the constrained optimization problem may be summarised as 
follows ; 

(1) The soluticffi is commenced from a feasible point 

A suitable value is chosen for the penalty parameter 
r and the counter k set equal to 1, 

(2) Using DFP algorithm and cubic interpolation method as 
explained in steps (6) through (11)/ the vector Xj^ is 
obtained which irdnimizes the modified objective func- 
tion |2f (X/ r^,). 


(3) The point Xj^ is tested for optimality. If it turns 
out to be optimal/ the process is terminated, 

(4) Otherwise/ the next value of peialty parameter is 
obtained frcxn the relation 


r, , = cr, / where c is less than 1. 

k+1 k 


( 5 ) The new value of k is set eq\aal to k+1 / the new 


C6) 




starting point X^ set equal to xj^ and the next mini- 
mization cycle starts from step (2), 

The DFP algorithm and the cubic interpolation 
method can be svimmarised as follows; 

The algorithm starts with the initial point X^ and a 
n X n positive definite symmetric matrix [h]^/ where 
n is the number of design variables. The iteration 
nuimber i is set equal to 1, 

Ihe gradient of the modified objective function/ 

0^/ at the point X^ is computed/ from which the 


( 7 ) 





search direction is found as 

= - [h]^ (6.1) 

iKr 

(8) Ihe rainiroizing step length ctj^ in the direction is 

foiind using the cubic interpolation method. The impro- 
ved design point is then four>d as 




Ar-. ' 

+ a,, s. 
1 1 


(6,2) 


(a) 


( 10 ) 


The new point is tested for optimality. If it is 
found to be optimal/ the iterative procedure is termi- 
nated , 


Otherwise/ the [h] . matrix is u:^ated as 

X 


[H]i^^ = [H]j^ + [m]^ + [n]^ 


-T 

e Si s_: 

where [M]. = ^ 

s: Q. 

11 


(6.3) 


(6.4) 


CnI 


([H]^ Qj.)([H3i Q^) 
^ [Hli 


(6.5) 


and 


^ i = ^J^i+1 - 


( 6 . 6 ) 


(11) The new iteration nvimber i is set equal to i+1 and the 
iteration started from st^ (7). 


6,3 Organisation of the Computer Programs 


The computer programs used in this study can be 
broadly classified into four categories; main program/ opti- 
mization subroutines/ function subroutines and the function 



subprograms. The main program essentially reads the input 
data, calls the relevant subroutines, and prints the output. 
The principal subroutines used to carry out the optimizaticn 
are MINIM4 and STEPX3, the former to find out the search 
directicai and the latter to obtain the minimizing step length. 
In order to find out the gradient of the modified objective 
function and the directional derivative, these subroutines 
in turn call GRADN and SL0PE2. The subroutine SLOPE 2 also 
calls GRADN. The subroutine GRADN makes use of the finite 
differoice approach to find the gradient vector. Since the 
gradiait, in the finite difference method, is computed on the 
basis of function value, subroutine GRADN in turn calls FIN. 
Ihe purpose of the subroutine PTN is to find the value of 
the objective fvuncticai, all the constraints g^ (j = 1^2/...m), 
and the modified objective function J2( for a chosai value of 
the penalty parameter. Thus, subroutine FIN in turn calls 
OBJT and CONST. The constraints can only be evalviated after 
carrying out the elastic and limit analysis of the tower. 

Thus, the subroutine CONST in turn calls ELANL and LIMANL. 
These in turn make use of the following subroutines: MATINV to 
obtain the inverse of a matrix for solving a set of linear 
simultaneous equations for finding the design forces in the 
reservoir portion, RAPSCX^ to solve a set of nonlinear simul- 
taneous equations for limit analysis of the cylindrical 
parti ticn, RECTI to find the ultimate strength of a rectan- 
gular section, and CIRCUl to obtain the ultimate strength of 
circular cross section of the supporting shaft. Function 
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subprograms are employed by subroutines RECTI and CIRCUl to 
find the value of functions for obtaining the values of 
neutral axis and mcroent of inertia, 

6.4 Some Important Aspects of Computer Programs 

Several precautions are necessary for obtaining the 
satisfactory restolts through canputer programs, particularly, 
vdien optimization technique is incorporated in the design 
process. Based upon the observations made during this study 
and the experience of others (Subramanyam, 1981; Srinivasa 
Rao, 1983), the following precautions are considered to be 
the important ones; 

(1) While carrying out the elastic analysis through force 
method, described in Section 3,3.1, it has been obser- 
ved that satisfactory results can only be obtained 

when double precision is employed during computations. 
When analysis is carried out adopting single precision 
it has been found that the values of stress resultants 
suddenly jump at some points, in some cases. The 
Table 6.1 shows the results of the analysis through 

force method for a tower with the following data; 

3 

capacity of the tank = 300 m , R 2 = 1.125 m, and 
R 1 /R 2 = 4,0. The table shows the values of stress 
resultants, along the wall of the conical tank, which 
are obtained using single as well a$ double precision. 
Thus, it can be seen that the solution obtained by 
single precision may lead to a serious error in the 
results. This happens due to the way in which the 
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Table 6.1 Results of analysis through force method using 
single and double precision 


T — — — 

; Values of stress resultants obtained using single 
precision (SP) and double precisicn (DP) 


u 

i in kN/m 

1 X 

i — 

[ Nq in kN/m 

in 

kN m/m 


: (SP) 

1 - 

i (DP) 

1 

i (SP) 

1 

j (DP) 

I (SP) 

1 - 

i (DP) 

D.OO 

-614.13 

-614.13 

163.16 

163.16 

1.9425 

1.9425 

D.05 

-511.33 

-511.33 

164.79 

164.79 

-0.083 

-0.083 

0,10 

-429.93 

-429.93 

184.49 

184.49 

-0.025 

-0.025 

0.15 

-363.88 

-363.88 

193.76 

193.76 

0.004 

0.004 

0.20 

-309.23 

-309.23 

201 . 24 

201.24 

-0.000 

-0,000 

0.25 

-263,22 

-263.22 

206.47 

206.47 

0.000 

-0,000 

0.30 

-223.99 

-223.99 

209.33 

209.33 

0.000 

0,000 

0.35 

-190.22 

, -190.22 

209.86 

209.86 

-0.000 

0.000 

o 

o 

-160.93 

-160.96 

208.04 

208.04 

0.026 

-0.000 

0,45 

-136.49 

-136.49 

209.75 

203.87 

0.000 

0.000 

0.50 

-113.27 

-113.27 

103.39 

197.37 

-0,000 

0,000 

0.55 

-105.11 

-93.88 

940.41 

188.52 

0.000 

0.000 

0.60 

-77.00 

-77.00 

3184.8 

177.34 

0.000 

0.000 

0.65 

-62.35 

-62.35 

163.81 

163.81 

439.75 

-0.000 

0,70 

-1955.3 

-49.73 

147.94 

147.94 

-0.000 

-0.000 

0.75 

-38.96 

-38,96 

129.72 

129.72 

0.000 

0.000 

0.80 

-29.90 

-29.90 

109.17 

109.17 

-0.000 

-0*000 

0.85 

-22.44 

-22.44 

86.25 

86.25 

-0.000 

-0.000 

0.90 

-16.48 

-16.48 

Very high 

61.06 

0.005 

0.005 

0,95 

-11.90 

-11.90 

34.42 

34.42 

-0.019 

-0.023 

1.00 

^8.75 

-8.75 

Very high 

-3.89 

-0.178 

-0.178 







Z.Z. / 


arithmetic operations are performed by the computer. 
Particularly in this case it is likely to happen because 
some of the elements of the flexibility matrix are of 
the order of 10 ® to Hence double precision 

should be employed vhile using force method of analysis 
for shell structures to overccane such difficulties. 

(2) All the ccaistraints should be normalized so that their 
values vary between ~1 and 0. If the constraints are 
not normalized then the contribution of the various 
constraints may be very much different in the formula- 
tion of the (ZT-f unctican . This will badly affect the 
convergence rate during the minimization of function. 
In view of the foregoing fact all the constraints in 
the present study have been expressed in their normal- 
ized form. 

(3) Too ro\ach difference in scale of the variables may cause 
some difficulties vhile selecting increments for step 
laigths or calculating nximerical derivatives. Some- 
times the objective function contours get distorted 
due to these scale disparities. Hence it is a good 
practice to scale the variables so that all the vari- 
ables become dimensionless and vary between a short 
range. However, the need of scaling the variables in 
the present study is not felt as the value of all these 
variables is expected to lie between 0.125 to 2,0 
approximately . It should be noted that the thicknesses 
of various elements are taken in metres in optimal 
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design process/ whereas^ these are given in mm in the 
design tables for conveni«ice. 

(4) The rate of convergence during the minimization of the 
0— function is found to be very much sensitive to the 
initial value of poialty parameter r^ and to some extent 
to reduction factor c. In the present study/ a value 

of r^/ which gives, the value of 0(Xj^, r^) approximately 
eqioal to 1.2 to 1.75 times the value of f(X^) is found 
to be quite satisfactory in most of the cases. Whereas / 
a value of 0.1 has been proihsd to be satisfactory for 
the reducticxi factor c in this study . 

(5) All the techniques available for the solution of non- 
linear programming problems guarantee/ at least, a 
local minimum. Ihis wotild ^Iso be the global mininHJm/ 
if the objective function is convex. It is not easy 
to ascertain whether the objective function is convex 
or not/ for the nature of functions encountered in this 
study or in other real life problems. To get over such 
a situation, the optimum solution is obtained with 
different starting points and the solutiais so ot^ained 
are compared. If the solutions match then it is very 
likely , but not certain , that the solutions correspond 
to the global minimum. The Table 6.2 shows the 
optimal designs of a water tower, obtained from three 
different starting points. It can be seen that there 
is no significant difference in the solutions obtained. 
In the present study, as all the problems solved are 



Table 6,2 Optimal design of a water tower from different starting points 

3 0 

Capacity = 500 m , “ 1*75 m, R-i/R^ = 4.0/ P^ = 1»5 kN/in^ 



230 


of the same natvire, the possibility of* having missed 
the global minimimi is practically nil. 

(6) In almost all constrained optimization problems^ there 
will be some constraints which are critically satisfied. 
Since the penalty term taids to infinity as the solu- 
tim approaches the constraint boundaries, special 
precautions have to be taken to overcome the overflow 
problem. 

(7) In order to avoid lannecessary large nxmber of vinconst- 
rained minimizaticnsas well as the prenature termina- 
tion, several convergence criteria have been employed 
in the optimization process. These include comparison 
of the relative values of the objective functions, the 
design variables, the gradient and related quantities, 
of the two successive iterations. In order to make use 
of the sol uti OTIS so obtained with confidence, and to 
test for optimality, perturbation of the design vector 
and testing the Kuhn-Tucker conditions can be resorted 
to. 



CHAPTER 7 


SUMMARY AND CONCLUSICNS 

7»1 Sunmairy 

Tbe main aim of this thesis has been to investigate 
the optimal limit state design of reinforced concrete, shaft 
supported, conical water tanks for Indian environmental condi- 
tions., The '«mjrtiasis has been to obtain practically usable 
optimal cost designs, conforming to relevant codes of practice, 
which can be useful to a professional engineer. 

The analysis and design have beai carried out using 
modem techniques which allow a useful insight into the 
Structxiral behaviour and design aspects of the structures 
c<^sidered in the present study under the expected loads. 
Pcongitiy and efficiency of these structures are ensured throuigh I 
optimization in which all the strength and serviceability I 

constraints of the relevant limit states are incorporated in | 
the formvLLated optimization problem. The optimization study 
of these water towers gives both optimal configuration and 
the optimal values of the design variables. Frarthermore, the 
influence of the escalation in unit costs of concrete, steel 
and formwork on the optimal configioration and the design 
variables of the tower has been studied in the present work. 

Parametric studies have been conducted for the most 
commcaily used reservoir capacities and staging heights for 
each of the three types of the water towers Shown in Figure 
1.1, The basic wind pressure is also treated as a parameter 



in the study of optimal designs of these water towers. The 
optimal designs so obtained are finally checked to ensure 
their safety and suitability for various seismic zones of 
India. 

Finally# the designs are presented in the tabular 
form, for eacdi of the three types of towers, fron which the 
optimal design variables and the parameters fixing the optimal 
configurati<»i can be picked up for a particular tank capacity, 
staging hei 5 ^ht, and wind and seismic zones. 

7,2 Conclusicns 

From the present study a number of conclusioas can 
be drawn. Ihese are listed as follows: 

(1) The use of limit state design philosophy and optimi- 
zaticn tedinique has resxalted in economical designs 
which remain safe and serviceable throughout their 
design life. 

(2) For every tank capacity, staging height, and wind and 

seianic zcaies there exist optimal values for the 
radius of the supporting shaft (^ 2 ^^ slope of 

the wall of the conical reservoir (a) with the hori- 
zontal for which the cost of the superstructure of 
the water tower is minim\jm. Similarly for water 
towers with given radius of supporting shaft, there 
is an optimal value of slope (“) which gives the 
minimum cost of the superstructure, 

(3) It is quite interesting to note that the optimal 
value of the slope (a) is found to lie within a 



very short range of 0.85 to 0,96, 0,991 to 1.02 and 
0.9 to 1,03 radians for Type I, Type II and Type III 
water towers respectively. Ihus, it can be concluded 
that the most probable optimal value of a will lie 
between 50“ to 60“ for all the cases. 

The optimal values of R 2 and a for a particular tank 
capacity, staging height, and wind and seismic zones 
in case of Type III towers are found to be less as 
compared to those of Type II water towers. In most 
of the cases. This is due to the fact that the cost 
of the partition in Type III towers reduces with 
decrease in the values of both R 2 and « . 

Ihe optimal configure tlcxi of the reservoir ad justs 
itself in a manner that the bending moments and shear 
forces at the junctions reduce considerably bringing 
the state of stress in the reservoir portion in the 
close proximity to that of the membrane state. It 
is also observed that the value of maximum hoop 
tension in the conical reservoir becomes the minimum 
for the optimal configuration, in most of the cases, 
Ihe location of the maximum hoop tension in the 
conical reservoir is invariably foiond to be between 
0.3 L to 0.5 L (Figure 3.5) from the junction 2, for 
all the cases of three types of towers considered 
in the present study. The variation of membrane 
stresses along the wall of the conical reservoir is 
also found to be similar in all the cases and the 



bending and shear forces remain limited to the boun- 
dary region of junctions only. 

For a specified cost ratio of concrete to steel, it 
is observed that the optimal, configuration and the 
design variables are not sensitive to the changes in 
the cost of formwork, whereas, for various cost 
ratios of ccancrete to steel, some of the optimal 
design variables do vary but slightly, while the 
optimal configuraticai remains essentially unaltered, 
Ibus for all practical purposes both optimal confi- 
gxuration and design variables remain unaffected by 
any change in the cost of materials. Hence the 
optimal configuration with corresponding optimal 
values of the design variables obtained by consid- 
ering the costs of materials prevailing at the time 
of investigation will remain optimal for all pract- 
ical purposes when the costs of materials are diff- 
erent , 

In Type III towers the provision of partition has 
resulted in considerable increase of banding moment 
values, for the various elements, in the boundary 
region of junction 2, This is observed, mostly, 
when only outer portion of the tank is full. How- 
ever, for the normal case, when both the portions 
are full, these bending moment values are comparable 
with those obtained for water towers without parti- 
tion. The governing design values of hoop tension 



for tile partition are obtained when inner portion of 
the tank is full , 

(9) The design forces in the reservoir portion are gover- 
ned by axisytnmetric loading combination of dead load 
and hydrostatic pressure loading under full tank 
ccaidition. Whereas, lateral forces due to wind and 
seismic effect under full or empty tank condition 
are found to be critical for the supporting shaft. 

(10) For towers with small capacity tank, the wind forces 
are found to predominate over the seismic forces for 
all staging heists ccxisidered in this study. It is 
observed that for 100 to 300 kl capacity tank a tower 
designed for wind zone I is found to be safe for all 
the seismic zones of India for almost all heights of 
staging ccaisidered herein. It is also identified 
that the ratio of seismic forces to wind forces 
increases with the increase in the capacity of the 
tank for a particular staging height and it decreases 
with the increase in the staging height for a givai 
capacity of the tank, 

(11) The serviceability limit state of cracking is the 
only one which is found to govern the design of the 
reservoir portion in all the cases with the excep- 
tion of some of the cases of 1000 ra^ capacity tank 
in which the constraint dealing with the xiltimate 
limit state of strength of a section of the conical 
tank at jmction 2 also becomes critical. In case 

3 

of supporting shaft, for towers with 100 and 200 ro 



capacity tank# only the serviceability limit state 
o£ deflection is found to govern the design. Whereas, 
for towers with larger capacity tank, both deflection 
as well as the liltimate strength of a local section 
of unit width of the supporting shaft become critical, 
Ihe resistance of these structures to collapse 
is invariably found to be more than adequate in all 
the cases, 

(12) Ihe extra cost for the provision of the partiticn 
for 500, 750 and 1000 kl capacity tanks has come out 
to be aiproximately 11 to 15%, 11 to 14% and 10 to 
12% respectively of the cost of the tower without 
partition. Ihe various benefits of partition, there- 
fore, accrue only by incurring the said extra expen - 
ditxire. 

(13) Fjrom the comparative cost study of two types of 
towers viz. Type I and Type II it is observed that 
for tank capacities upto 300 kl the difference 
between the costs of the two towers is marginal (1 to 
3,5%). However, for capacities 500 kl and above 

the Type II towers are quite economical. Although 
the cost of the Type I towers for capacities upto 
300 kl is slightly higher than the Type II towers, 
from constructional point of view the former are 
simpler. Hence, in the present study, the prefer- 
ence of Type I towers over Type II upto 300 kl 
capacity is justifiable. 
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(14) As the present study addresses itself to the super- 
striKiture only/ the optimal designs for radii of 
supporting shaft other than the optimal radius are 
also givQi in tabular form. These designs are in 
additicxi to that corresponding to the optimal radius 
ai^ may come out to be economical for a particular 
foundation condition vdien the cost of the foundation 
is also takai into consideration. 

(15) Althoxogh the present study is confined to Indian 
aivircaunental caiditicxis^ the results can be judici- 
ously utilized for other environments as well. 

7.3 Sxiggestions for Future Work 

The following are some of the suggestions which may 

be considered for future work: 

(1) Study of optimal designs of overhead tanks of different 
shapes for different capacities and staging heists 
will help in choosing the most economical shape, 

(2) The extension of present work for optimal design of 
overhead conical water tanks with a partition providing 
a definite capacity of the inner portion is straight- 
forward, 

(3) Although# for larger capacity tanks# an indirect 
procedure given in the present study for tapering the 
thickness of the wall of conical reservoir works 
quite satisfactorily and results in considerable 
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saving of volume of concrete, the optimal designs with 
tapering thickness of the reservoir using an exact 
method of analysis can be investigated. 
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